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ABSTRACT 


The  physical  characteristics  of  conifers  that  had  been 
fed  upon  by  porcupines  (Erethizon  dors a turn)  and  of  unused 
conifers  were  examined  on  two  areas  of  montane  forest  in 
southwestern  Alberta.  The  three  vegetation  types  examined 
included  a  windward  community  co-dominated  by  limber  pine 
and  Douglas  fir;  a  leeward  forest  of  Douglas  fir  with 
scattered  hybrid  spruce;  and  a  leeward  community,  also  co¬ 
dominated  by  limber  pine  and  Douglas  fir,  which  appeared  to 
be  serai  to  the  Douglas  fir  forest.  The  average  age  of 
conifers  in  these  forests  was  approximately  65  years. 

Of  the  factors  considered,  exposure  and  age  were  most 
influential  in  determining  the  size  and  form  of  trees 
measured  on  both  areas,  although  density  was  also  important 
in  leeward  forests.  The  superior  growth  rates  of  leeward 
trees  relative  to  windward  trees  and  the  associated 
differences  in  phloem  production  probably  account  for  the 
preference  of  porcupines  for  leeward  trees.  Although  spruce 
grew  more  rapidly  than  Douglas  fir  and  individuals  of  both 
species  grew  faster  than  windward  limber  pine,  porcupines 
did  not  exhibit  apparent  interspecific  preferences. 

The  trees  used  by  porcupines  commonly  had  extensive 
crowns  and  were  larger  and  more  vigorous  than  unused  trees, 
so  that  they  probably  provided  greater  access  to  more 
abundant,  better  quality  food.  This  type  of  tree  was  most 
frequent  in  low  density  stands  of  the  leeward  forests. 
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A  greater  proportion  of  the  conifers  on  the  Zoratti 
area  were  used  by  porcupines*  even  though  trees  on  the 
Skelton  area  were  more  vigorous.  This  difference  between 
areas  appears  to  have  been  associated  with  the  proximity  of 
the  Zoratti  forests  to  summer  feeding  areas  rather  than 
differences  in  food  or  habitat  quality.  Because  of  the 
large  number  of  porcupines  using  the  Zoratti  area*  the 
average  size  of  used  trees  was  smaller  than  on  the  Skelton 
area.  In  addition*  the  pattern  of  repeated  use  of  Douglas 
fir  in  the  leeward  pure  forests  differed  between  areas. 

The  use  of  both  areas  by  porcupines  peaked  and  declined 
in  relation  to  the  ages  of  the  forests*  suggesting  that 
trees  may  become  too  large  to  be  climbed  and/or  the 
declining  growth  rates  of  trees  are  associated  with  a 
reduction  in  the  quality  of  inner  bark.  During  winters  of 
porcupine  abundance*  the  use  of  poor  quality  stands 
increased.  The  nature  of  the  observed  preferences  indicates 
some  active  selection  of  trees  by  porcupines.  Although  the 
actual  mechanisms  of  tree  selection  are  unknown*  stand 
density  and  the  physical  characteristics  of  individual 
trees*  particularly  trunk  diameter  and  crown  size*  are 
correlated  with  porcupine  preferences.  Winter  food  and 
habitat  selection  by  porcupines  in  these  areas  appear  to  be 
interdependent. 
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INTRODUCTION 

Environments  are  not  homogeneous.  If  the  relative 
ability  of  a  habitat  to  provide  food  and  protection,  during 
both  reproductive  and  fallow  periods,  is  taken  as  a  measure 
of  habitat  quality,  those  animals  that  can  consistently 
detect  and  utilize  the  higher  quality  portions  of  their 
environment  should  make  the  greatest  contribution  to  future 
generations.  Variation  in  the  value  of  food  (net 
metabolizable  nutrients)  available  within  a  particular 
habitat  should  increase  the  scope  of  selectivity  that  could 
be  exercised  by  a  given  consumer  to  optimize  fitness 
(Eadie  197  0).  Food  selection  should  therefore  be  an 
integral  part  of  the  habitat  utilization  pattern  of  any 
mobile  animal. 

Porcupines  ( Erethizon  dor sat um)  are  particularly  good 
subjects  for  examining  the  association  between  food  and 
habitat  selection  because  evidence  of  their  winter  feeding 
is  both  obvious  and  long-lived.  During  the  summer, 
porcupines  feed  on  herbaceous  vegetation  and  the  foliage  and 
catkins  of  some  deciduous  trees  (Dodge  1967;  Brander  1973). 
Concurrent  with  the  cessation  of  plant  growth  and  the  onset 
of  breeding  in  the  fall,  porcupines  move  to  areas  more 
heavily  wooded  than  their  summer  ranges  and  begin  to  feed  on 
the  inner  bark  of  woody  plants  and  the  foliage  of  several 
coniferous  species  (Gabrielson  and  Horn  1930;  Curtis  and 
Wilson  1953;  Brander  1973).  When  feeding  on  inner  bark, 
porcupines  remove  both  the  cambium  and  the  phloem,  exposing 
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the  xylem.  Daring  subsequent  growth  periods  this  exposed 
xylem  is  slowly  covered  by  the  centripetal  growth  of  callus 
tissue  formed  from  the  surrounding  cambium  (Brown  1971b). 

The  resulting  scars  are  a  semi-permanent  record  of  the 
amount  and  location  of  bark  removed  and  of  the  year  in  which 
feeding  occurred.  Although  barking  may  occur  at  ground 
level,  especially  in  saplings,  most  feeding  occurs  within 
the  crown  of  the  tree  (Taylor  1935)  •  Since  porcupines  often 
use  the  same  tree  for  food  and  shelter,  winter  food  and 
habitat  selection  by  this  species  may  be  intimately  related. 

The  preferential  use  of  a  portion  of  a  resource  need 
not  be  a  direct  result  of  selective  behaviour  because  the 
dispersion  of  a  group  of  consumers  may  be  controlled  by  the 
distribution  of  other  essential  resources  (Wiens  1976) • 
Porcupines  could  exhibit  preferences  at  four  successively 
finer  subdivisions  of  a  forest:  forest  types;  stands  within 
a  type;  individual  trees;  and  portions  of  a  free.  The 
quality  of  each  of  these  subdivisions  is  dependent  upon  the 
quality  of  every  other  subdivision  as  well  as  the  direct 
influences  of  the  physical  environment  (Billings  1952) ,  and 
it  was  the  effect  of  these  interrelationships  on  porcupine 
preferences  which  formed  the  focus  of  this  study.  Of 
particular  interest  were  characteristics  of  feeding  that 
were  suggestive  of  active  food  or  habitat  selection. 

Previous  ecological  studies  of  porcupines  have  been 
largely  directed  at  the  silvicultural  significance  of  this 
species  (Taylor  1935;  Curtis  1941;  Krefting  et  al.  1962; 
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van  Deusen  and  Myers  1962)  .  Although  these  studies  are  the 
source  of  much  of  the  information  regarding  porcupines, 
their  ultimate  interests  were  with  the  forest  rather  than 
with  the  animals.  An  alternative  approach,  which  emphasized 
the  dependence  of  porcupines  on  the  nature  of  the  forest, 
was  adopted  during  this  study. 
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STUDY  AREAS 

Two  areas  were  intensively  examined  during  the  course 
of  this  study,  both  located  on  the  eastern  fringes  of  the 
foothills  belt  that  runs  along  the  base  of  the  Front  Ranges 
of  the  Rocky  Mountains,  The  Zoratti  study  area 
(49°2  7  *  N;  1 1 4°0  4 '  W :  Sec,  10,  15,  16  Tp.  6  R.  1  W,  5), 
located  11  km  west  of  Pincher  Creek,  Alberta,  covers  an  area 
of  36.3  ha  and  varies  in  elevation  from  1210  to  1320  m 
(Figure  1).  The  Skelton  study  area  (49°349  N;  114°139  W: 
Sec.  21,  22  Tp.  7  R.  2  W.  5) ,  situated  16  km  northwest  of 
the  Zoratti  site,  encompasses  115.4  ha  and  ranges  in 
elevation  from  1180  to  1360  m  (Figure  2). 

Parallel  ridges  of  thrust-faulted  and  folded  Cretaceous 
sandstones  with  a  northwest-southeast  orientation 
characterize  the  foothills  of  this  region.  Fault  blocks  in 
the  area  were  thrust  over  adjacent  blocks  to  the  east  so 
that  west-facing  slopes  are  gentler  than  those  facing  east. 
With  the  exception  of  portions  of  the  ridgetops,  the  entire 
region  is  covered  by  a  veneer  of  glacial  till  of  varying 
thickness  (Beaty  1975) .  Although  both  study  areas  were 
located  on  this  system  of  ridges,  exposure  on  the  Skelton 
site  was  dominated  by  the  valley  of  the  Crowsnest  River  to 
the  north.  This  valley  was  carved  by  an  early-Pleistocene 
valley  glacier  which  cut  across  the  system  of  ridges, 
greatly  altering  local  aspect  (Beaty  1975)  . 
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Figure  1. 


Map  of  the  Zoratti  study  area  and  its  geographical 
location  (inset) .  Contour  interval  10  m. 
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Figure  2.  Map  of  the  Skelton  study  area  including  the 

distribution  of  vegetation  types  and  the 
locations  of  the  sampling  plots  (open  circles) . 
Contour  interval  10  m. 
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The  proximity  of  this  region  to  the  Crowsnest  Pass  has 
a  marked  effect  on  its  climate.  Opening  directly  onto  the 
plains,  this  pass  serves  as  a  channel  for  Pacific  air 
flowing  into  Alberta  (Longley  1972).  As  a  result, 
prevailing  westerly  winds  are  often  intense  and  seemingly 
omnipresent  (Table  1).  These  winds  often  take  the  form  of  a 
Chinook  (foehn) ,  involving  a  warm,  dry  air  mass,  which 
greatly  increases  evapotranspiration  in  exposed  locations. 
The  Pincher  Creek  region  can  expect  almost  300  cm  of  snow 
during  the  winter  (Table  2) ,  but  because  of  the  frequency  of 
warm,  westerly  winds  the  exposed  slopes  of  the  foothills  are 
bare  much  of  the  winter,  while  the  sheltered  leeward  slopes 
are  snow-covered. 

The  physiognomy  and  climate  of  the  foothills  west  of 
Pincher  Creek  are  together  responsible  for  a  distinct 
distribution  of  plant  communities  (Figure  3:  Moss  1944).  A 
brief  description  of  the  four  major  vegetation  types 
recognized  on  the  two  study  areas  (Figures  2  and  4)  follows; 
more  detail  can  be  found  in  Moss  (1944)  and  Moss  and 
Campbell  (1947) . 

A  Festuca-Danthonia  community  (Moss  1944)  ,  regarded  by 
Moss  and  Campbell  (1947)  as  an  edaphic  climax  within  the 
Festuca  scabrella  association  maintained  by  moderate 
grazing,  dominated  south-  and  west-facing  slopes. 

An  open,  limber  pine  (Pin us  f lexilis) -Douglas  fir 
(P seudo t suga  menziesii  var.  glauca)  woodland  occupied  the 
ridgetops  and  the  more  windswept  western  exposures.  The 
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Table  1.  Percentage  frequency  and  mean  speed  of  winds  from 
the  most  prevalent  directions  at  Pine her  Creek, 
Alberta,  1960-1972,  information  supplied  by  the 
Atmospheric  Environment  Service  of  Environment 
Canada. ) 

Direction 

SW  WSH  W  WNW  Calm 


Month 

% 

km/hr 

% 

km/hr 

% 

km/hr 

% 

km/hr 

% 

January 

4 

19.5 

15 

36,5 

37 

35.5 

6 

20.0 

15 

February 

4 

18.4 

18 

34.  6 

40 

33.  0 

5 

18.  1 

14 

March 

6 

18.7 

20 

29.9 

33 

28.3 

5 

16.0 

14 

April 

7 

19.7 

23 

30,4 

30 

27.7 

5 

14.9 

12 

May 

7 

18.2 

13 

26.  1 

30 

24.2 

6 

14.9 

12 

June 

7 

16.6 

22 

26.2 

29 

22.  7 

7 

16.2 

14 

July 

9 

16.6 

18 

22.6 

28 

21.1 

7 

15.0 

16 

August 

9 

17.4 

18 

23.5 

31 

21.3 

7 

14.  1 

15 

September 

8 

17.3 

20 

25.  6 

31 

23.8 

6 

15.  4 

15 

October 

6 

19.0 

24 

31.4 

41 

30.  1 

6 

18.2 

10 

November 

5 

19.0 

21 

31.4 

40 

31.4 

6 

19.8 

13 

December 

4 

20.0 

16 

33.8 

41 

32.6 

7 

20.5 

15 

' . 
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Table  2.  Temperature  and  precipitation  means  for  the 

winter  months  at  Pincher  Creek,  Alberta,  1960-1972. 
(Information  supplied  by  the  Atmospheric 
Environment  Service  of  Environment  Canada.) 

Temperature  (°C)  Precipitation  (cm) 


Month 

Daily 

Maximum 

Daily 

Minimum 

Daily 

Mean 

Total 

Snow 

October 

12.  9 

0.6 

6.8 

2.79 

25.  1 

November 

4.8 

-6.3 

-0.7 

2.59 

30.0 

December 

-1.6 

-12.2 

-6.8 

3.33 

42.9 

January 

-3.9 

-15.5 

-9.7 

3.89 

52.3 

February 

2.2 

-9.3 

-3.6 

2.44 

32.0 

March 

3.  9 

-8.1 

-2.1 

3.28 

41.9 

April 

9.8 

-2.  4 

3.7 

6.12 

66.5 

\ 
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Figure  3.  A  profile  of  the  Zoratti  study  area,  taken  along 

the  transect  indicated  by  the  opposing  pointers 
in  Figure  4,  illustrating  the  relationship 
between  slope  and  exposure  and  the  distribution 
of  vegetation. 
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Figure  4,  flap  of  the  Zoratti  study  area  including  the 

distribution  of  vegetation  types  and  the 
locations  of  the  sampling  plots  (open  circles) . 
The  left-right  axis  has  been  expanded  by  a  factor 
of  five.  Contour  interval  10  m. 
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relative  proportions  of  the  two  tree  species  varied  with  the 
degree  of  exposure  and  soil  depth,  limber  pine  being  more 
common  on  the  least  sheltered  sites  with  exposed  bedrock. 
Most  trees  in  this  community  exhibited  the  stunted, 
asymmetrical  form  typical  of  trees  growing  near  timberline. 

Douglas  fir  forest  occupied  the  leeward  slopes  of  both 
study  areas.  Hybrid  spruce  fPicea  glauca  x  P.  engelmannii ; 
Daubenmire  1974)  are  a  minor  component  of  this  community, 
and  other  coniferous  species,  such  as  lodgepole  pine 
(Pin us  con tort a  var.  latif olia)  and  subalp in e  fir 
(Abies  lasiocarpa) ,  are  rare.  Howe  (1972)  stated  that  the 
montane  forest,  of  which  this  community  is  a  part,  developed 
in  response  to  the  prevailing  dry  climate,  although  the 
lower,  and  presumably  eastern  distributional  limit  of 
Douglas  fir  is  influenced  by  summer  moisture  deficiencies 
(Daubenmire  1943) .  Areas  to  the  west  of  the  study  sites 
receive  more  precipitation  (Wyatt  et  al.  1939)  and  support 
extensive  stands  of  hybrid  spruce  characteristic  of  the 
subalpine  forest  of  the  eastern  slopes  of  the  northern  Rocky 
Mountains  (Rowe  1972). 

With  decreasing  slope  along  the  base  of  the  ridge,  the 
leeward  Douglas  fir  forest  on  the  Zoratti  study  area  graded 
into  a  younger  limber  pine-Douglas  fir  community  that  is 
apparently  serai  to  a  purer  Douglas  fir  forest.  The 
presence  of  limber  pine  as  a  pioneering  tree  species  may 
have  been  a  result  of  the  isolation  of  this  area  from  the 
western  subalpine  forest.  The  usual  pioneer  of  eastern 
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slope  forests  is  lodgepole  pine  (Cormack  1953) ,  but  only  a 
single  lodgepole  pine  was  found  on  this  study  area. 

Aspen  poplar  (Pod ulus  tremuloides)  groves  are  located 
in  flatter  areas  which  are  well  drained  and  protected  from 
the  wind.  As  the  quality  of  drainage  declines,  the  aspen 
give  way  to  balsam  poplar  (P.  balsam if era)  and  then  to 
willow  tSalix  spp.)  and  water  birch  (Betula  occidentalis) . 
Willows  {5«  bebbiana,  S.  scouleriana,  and  others)  are  the 
only  local  deciduous  species  that  are  used  as  winter  food  by 
porcupines. 

Although  a  few  fire-scarred  trees  exceeding  250  years 
of  age  are  known  from  both  study  areas,  the  vast  majority 
are  less  than  90  years  old.  Dawson  (1886)  reported  large 
tracts  of  burned  timber  in  the  foothills  of  southern  Alberta 
which  he  attributed  to  the  initial  use  of  the  Crowsnest  Pass 
by  white  travellers  during  the  latter  half  of  the  nineteenth 
century.  The  forests  on  both  study  areas  are  probably  the 
products  of  regeneration  following  extensive  fires  at  that 


time. 


'  •:  "<  *•  -  ■*  •  v  ■  3 

,  _ .  •  .  v  ..  *  > 

■ 

. 


18 


METHODS 


A  total  of  2561  trees  were  sampled  on  20  x  20  m  plots 
during  the  summers  of  1974  (Zoratti)  and  1975  (Skelton). 

The  38  plots  examined  on  the  Zoratti  area  (Figure  4)  were 
regularly  distributed  along  transects,  separated  by  100  m 
intervals,  that  traversed  the  ridge  perpendicular  to  its 
long  axis.  Centers  of  adjacent  plots  on  the  same  transect 
were  45  m  apart  so  that  10  percent  of  the  total  area  of  the 
coniferous  forests  was  sampled.  A  preliminary  analysis  of 
the  data  collected  from  these  plots  suggested  that  the 
intensity  of  porcupine  use  might  have  been  related  to  stand 
density.  To  examine  this  relationship  more  rigorously,  the 
40  plots  sampled  on  the  Skelton  area  were  chosen  from 
190  plots,  distributed  as  on  the  Zoratti  area,  on  the  basis 
of  the  number  of  trees  per  plot  and/or  the  uniformity  of 
tree  distribution  (Figure  2) .  Several  plots  encompassing 
fewer  than  20  trees  were  enlarged  to  cover  an  area  of 
500  sg.  m.  Stand  density  (number  of  trees/100  sq.  m)  was 
calculated  as  though  each  plot  was  level.  Slope  and 
exposure  were  measured  for  all  plots  using  a  compass  and  a 
Suunto  clinometer. 

Species,  circumference  at  breast  height  (±  0.6  cm), 
height  (±  0.3  m  for  trees  under  7  m;  ±  1.5  m  for  taller 
trees) ,  and  evidence  and  the  probable  cause  of  leader  damage 
were  determined  fcr  all  trees  taller  than  1.35  m  (breast 
height) .  A  random  selection  of  20  percent  of  the  trees 
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within  each  plot  were  aged  using  an  increment  borer  and  a 
IQx  hand  lens.  Cores  were  extracted  at  stump  height  since 
the  age  differential  between  conspecific  individuals  in 
different  communities  would  be  less  at  that  height  than  at 
breast  height.  Each  tree  was  then  climbed  and  inspected  for 
porcupine  feeding  scars,  and  a  record  made  of  the  intensity 
of  use  (four  classes  based  on  the  total  scar  area)  for  the 
tree  as  a  whole  and  separately  for  the  lateral  branches  and 
the  bole.  In  addition,  the  age  of  use  was  established  and 
the  height  range  and  intensity  of  use  for  each  age 
estimated.  Ose  age  was  ascertained  by  removing  a  portion  of 
the  scar  tissue  surrounding  the  feeding  area  with  a  knife 
and  counting  the  annual  rings  from  the  level  of  exposure  by 
the  porcupine  to  the  most  recent  xylem  surface.  This 
technique  was  a  modification  of  the  procedure  employed  by 
Spencer  (1964)  to  document  porcupine  population  fluctuations 
over  a  225  year  period. 

Although  Douglas  fir  and  limber  pine  were  sufficiently 
abundant  to  be  adequately  sampled  with  plots,  spruce  were 
poorly  represented  on  both  study  areas  and  were  generally 
found  in  clumps  within  the  pure  fir  forest.  Since  the 
relative  preferences  for  these  species  were  of  interest,  the 
leeward  pure  forest  on  the  Zoratti  area  was  systematically 
searched  for  spruce  trees.  One  tree  was  randomly  selected 
from  every  three  encountered  and  examined  in  the  same 
fashion  as  trees  within  the  plots. 
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Several  parameters,  not  measured  directly,  were 
calculated  for  each  tree.  Diameter  at  breast  height 
(DBH  ±  0.2  cm)  was  derived  from  the  measured  circumference. 
The  ratio  of  the  DBH  to  the  height  served  as  an  index  of 
growth  form,  based  on  the  direct  relationship  between  bole 
diameter  and  crown  width  (Mitchell  1969) .  This  conformation 
index  (Cl),  generally  ranging  from  0.75  to  3.00,  was  only 
useful  for  intraspecific  comparisons  because  growth  form 
varies  between  species. 

Upon  examination  of  the  data  collected  during  1974,  it 
became  apparent  that  the  measured  diameter  at  breast  height 
would  only  be  indirectly  related  to  the  dimensions 
considered  by  porcupines  during  previous  years.  To 
circumvent  this  problem  an  attempt  was  made  to  determine  the 
DBH  and  height  of  each  tree  for  the  preceding  20  years. 

Cores  extracted  during  the  ageing  of  trees  on  the  Skelton 
study  area  were  placed  in  plastic  milkshake  straws  until 
they  could  be  inspected  in  the  laboratory.  In  addition,  the 
thickness  of  the  bark  (±  0.025  mm),  extracted  at  breast 
height  with  the  increment  borer,  was  measured  using  a  pair 
of  vernier  calipers.  Since  the  cores  extracted  on  the 
Zoratti  study  area  had  been  discarded,  this  area  was 
revisited  and  cores  collected  from  165  trees  randomly 
selected  from  the  sample  plots. 

In  the  laboratory,  the  annual  growth  increment 
l±  5  x  10~*  mm)  for  the  past  20  years  of  each  core  was 
measured  with  the  aid  of  a  micrometer  eyepiece  and  a  Wild 
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dissecting  microscope  (200x)  .  This  process  was  aided  by 
cutting  a  flat  surface,  perpendicular  to  the  long  axis  of 
the  tracheids,  with  a  scalpel  and  viewing  the  cores  under 
reflected  light.  The  length  of  each  core,  from  the  pith  to 
its  outer  surface,  was  measured  with  vernier  calipers. 

Calculation  of  the  DBH  of  any  cored  tree,  x  years  prior 
to  measurement,  employed  the  following  basic  algorithm.  The 
radius  inside  the  bark  during  the  year  of  interest  (IRx)  was 
taken  as  the  product  of  the  inside  radius  at  the  time  of 
measurement  ( (DBH/2) -BARK)  and  the  proportion  of  total 
length  (TC)  represented  by  the  length  of  the  growth 
increment  x  years  previous  (TC-25  In,  where  In  is  the  length 
of  the  growth  increment  for  any  given  year,  n) .  Therefore, 

I Rx=[  (DBH/2) -BARK  ]•  l  (TC-5  In)/TC)  (1) 

and  the  DBH  during  that  year  was  calculated  using  the  linear 
regression  equation  relating  the  inside  radius  at  the  time 
of  measurement  to  the  measured  DBH  (Husch  et  al.  1972) .  The 
error  associated  with  this  regression  equation  is  probably 
negligible  since  the  smallest  coefficient  of 
determination  (r2)  calculated  was  0.997. 

Since  cores  were  extracted  at  stump  height,  a 
correction  term  was  required  to  account  for  the  shorter 
total  core  length  at  breast  height.  The  age  of  a  tree  at 
breast  height  (ABH)  was  estimated  using  regression 
equations,  of  the  form  log  (HT) =a+b  (1/AGE)  (Husch  1963), 
relating  height  and  age  for  each  species  in  each  community. 
The  length  of  the  proximal  portion  of  the  increment  cores 
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representing  growth  prior  to  this  estimated  age  was  then 
subtracted  from  the  total  core  length  providing  an  estimate 
of  the  appropriate  core  length  (TC)  for  calculating  past 
diameters.  An  examination  of  the  error  associated  with  TC 
is  presented  in  Appendix  1. 

Using  these  data  for  cored  trees  from  each  forest  type, 
multiple  regression  equations  describing  the  relationship 
between  the  calculated  DBH  and  the  age  and  measured  DBH  were 
computed  for  each  of  the  20  years  preceding  measurement. 
These  equations  were  subsequently  used  to  estimate  the  past 
diameters  of  uncored  trees.  Age  was  included  in  all 
20  prediction  equations  only  if  it  was  a  significant 
contributor  to  at  least  half  of  the  series;  otherwise  the 
measured  DBH  was  used  as  the  sole  predictor  in  all 
equations.  The  use  of  age  in  these  equations  required  an 
estimate  of  age  for  each  tree  derived  from  the  appropriate 
linear  regression  equation  relating  the  age  and  DBH  of  cored 
trees. 

The  procedure  applied  in  estimating  the  height  of  any 
tree,  during  the  20  years  prior  to  measurement,  utilized  the 
relationship  between  height  and  DBH.  Of  the  four  general 
equations  suggested  by  Husch  et  al.  (1972) ,  the  power  curve 
(log-log  transformation)  provided  the  best  fit  to  the 
regression  of  the  measured  height  on  the  measured  DBH.  The 
growth  of  any  tree  was  therefore  taken  as  a  straight  line, 
on  a  log-log  plot,  between  the  measured  height  and  DBH  and  a 
height  of  1.35  m  and  zero  DBH,  allowing  the  estimation  of 
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the  height  corresponding  to  any  previously  established  DBH. 
Each  value  in  this  calculation  was  increased  by  one  to  avoid 
the  technical  problem  involved  in  determining  the  log  of 
zero  (Sokal  and  Rohlf  1969),  If  the  height  growth  had  been 
interrupted  by  a  factor  other  than  a  porcupine  killing  the 
leader  or  if  such  porcupine  damage  occurred  greater  than 
twenty  years  before  measurement,  previous  heights  were  not 
calculated*  Had  the  leader  been  killed  by  porcupines  within 
20  years  of  the  time  of  measurement,  the  growth  since  the 
use  occurred  was  estimated  using  a  regression  line  between 
the  height  of  the  top  of  the  uppermost  feeding  scar  and  the 
calculated  DBH  for  the  year  during  which  feeding  occurred 
and  the  measured  height  and  DBH* 

Although  errors  in  measuring  tree  and  use  age  and 
yearly  growth  increment  could  have  arisen  by  overlooking 
false  and/or  missing  rings  (McGinnies  1963) ,  this  type  of 
error  was  probably  insignificant.  Douglas  fir  and  limber 
pine  are  both  considered  excellent  dendrochronological 
material,  in  part  because  of  a  high  correlation  between  ring 
formation  and  the  annual  climatic  cycle  (Schulman  1956) . 
Schulman  examined  Douglas  fir  from  three  southwestern 
Alberta  sites  and  found  the  frequency  of  locally  absent  bole 
rings  ranged  from  0*1.1  percent  of  the  total  number  of 
rings,  depending  upon  the  degree  of  moisture  stress.  He 
also  stated  that  ”...  false  annual  rings  are  of  rare 
occurence  in  the  dendrochronologic  species  of  the  northern 
Rocky  Mountains....”  (pg.  19).  The  age  of  very  old  scars 
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located  on  senescent  branches  at  the  base  of  the  crown  was 
probably  less  accurately  measured  than  tree  age  because 
these  old  branches  often  expend  more  energy  in  respiration 
than  they  fix  through  photosynthesis  and  xylem  rings  are 
produced  irregularly  (Reukema  1959) . 

A  variety  of  statistical  techniques  was  applied  in 
describing  and  interpreting  the  relationships  between 
measured  variables  and  in  data  reduction.  Most  of  these 
tests  are  in  common  usage  and  do  not  warrant  detailed 
explanation.  The  multivariate  techniques  are,  however,  not 
as  widely  known  and  a  brief  description  of  their  major 
features  has  been  provided  in  Appendix  2.  Discussion  of  the 
statistical  procedures  followed  in  examining  specific 
problems  generally  precedes  the  presentation  of  the  related 
results.  The  log  likelihood  ratio  test  (G)  was  used  in 
preference  to  the  chi-square  test  for  the  reasons  outlined 
by  Sokal  and  Rohlf  (1969)  .  A  probability  level  of  <  0.05 
was  considered  statistically  significant  for  all  tests. 

Different  parametric  tests  were  used  in  the  comparison 
of  two  means,  depending  upon  the  relative  sizes  of  their 
respective  variances  as  determined  by  an  F-test  (Sokal  and 
Rohlf  1969).  The  statistic  derived  from  two  samples  with 
homogeneous  variances  is  based  on  a  pooled  variance  estimate 
and  is  distributed  as  Student's  t  with  (nl  +  n2  -  2)  degrees 
of  freedom.  If  the  variances  were  heterogeneous  the  pooled 
variance  estimate  becomes  meaningless  and  the  test  statistic 
is  therefore  calculated  using  separate  variance  estimates 
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and  is  approximately  distributed  as  Student's  t  with 

_ [  (Sl2/n1)  +  ( S22/n2 )  j2 _  ;  rounded  to  the 

[  (Sl2/n1)  2/ (nl-1)  ]+[  (S22/n2)  2/ (n2-1)  ]  nearest  intege 

degrees  of  freedom  (Nie  et  al.  1975) .  As  a  result  of  these 

differences  all  references  to  these  tests  will  include:  the 

test  statistic  (t) ;  the  probability  of  a  similar  random 

difference;  the  type  of  variance  estimate  used  (S=separate; 

P=pooled) ;  and  the  appropriate  degrees  of  freedom. 

Calculating  past  tree  dimensions  and  the  subsequent 
statistical  analysis  would  not  have  been  possible  without 
the  facilities  of  the  Department  of  Computing  Services  of 
the  University  of  Alberta.  The  Statistical  Package  for  the 
Social  Sciences  (Nie  et  al.  1975),  supplemented  by  several 
APL  programs,  provided  all  of  the  statistical  tests 
required.  Problem-specific  Fortran  programs  were  used  for 
data  manipulation. 
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RESULTS  AND  DISCUSSION 

Characteristics  of  Available  Trees 

If  the  physical  size  of  a  tree  influenced  its  selection 
as  a  feeding  tree  by  a  porcupine,  the  attractiveness  of  a 
particular  forest  would  be  in  part  determined  by  the  number 
of  available  trees  with  the  preferred  proportions.  The 
appeal  of  this  forest  would  be  further  influenced  by  the 
relative  quality  of  adjacent  forests,  so  that  the  factors 
affecting  tree  growth  and  distribution  would  indirectly 
affect  the  selection  of  food  and  habitat  by  porcupines. 

Since  a  tree's  dimensions  change  with  age,  the  appeal  of  a 
forest  or  tree  could  also  be  expected  to  change  through 
time. 

Although  exposure  was  varied  on  the  Skelton  area,  the 
open  limber  pine-Douglas  fir  forest  (windward  community) 
occurred  predominantly  on  west-facing  slopes  while  the 
closed  canopy  Douglas  fir  forest  (leeward  community) 
occupied  more  northern  exposures  (Table  3) •  Since  the 
measured  exposures  for  these  two  communities  were 
symmetrically  distributed  around  their  respective  means 
(mean,  mode,  and  median  vary  by  only  three  degrees  for  both 
communities) ,  these  mean  exposures  have  been  considered  the 
optimum  positions  for  each  vegetation  type.  Only  two 
exposures  existed  on  the  Zoratti  area,  neither  of  which 
approached  these  suggested  optimum  positions,  so  that  the 
equivalent  leeward  forest  (leeward  pure)  on  the  Zoratti  area 
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Exposure,  slope,  age,  and  species  composition  of 
the  conifer  dominated  communities  on  the  Zoratti 
and  Skelton  study  areas.  Each  entry  includes  the 
mean,  standard  deviation  and  number  of  observations. 


Exposure  Species 

(degrees  and  Slope  Age  Composition 

direction)  (degrees)  (years)  (%  Douglas  fir) 


Windward 


Zoratti  (X) 

220.0  SW 

12.5 

67.0 

16.5 

(S) 

— 

6.93 

26.71 

18.86 

(n) 

18 

18 

137 

1  1 

Skelton 

265.0  W 

15.8 

74.4 

68.9 

35.00 

9.30 

19.31 

34.79 

11 

11 

76 

1  1 

Leeward  Pure 

Zoratti 

40.0  NE 

34.6 

60.5 

97.9 

5.59 

19.31 

4.10 

13 

13 

112 

13 

Skelton 

357.1  N 

20.1 

66.5 

98.6 

33.70 

6.37 

16.08 

5.98 

29 

29 

246 

29 

Leeward  Mixed 

Zoratti 

40.0  NE 

24.9 

49.0 

49.4 

— 

5.08 

18.77 

26.80 

6 

6 

55 

6 

Combined 

Zoratti  61,6 

23.70 

304 

68.4 
17.19 
322 


Table  3. 


Forest 

Type 


Skelton 


* 

• 

• 

« 

• 

, 

0 

. 

* 


28 


was  restricted  to  the  steeper  slopes  (t=7.044  p<0.001  P 
df=40)  which  receive  less  solar  radiation.  The  younger 
(t  =  3.657  p<0.001  P  df=165)  mixed  forest  (leeward  mixed) 
described  earlier,  is  colonizing  the  gentler  slopes  (t=3.575 
0.01>p>0.001  P  df=16)  on  the  leeside  of  the  Zoratti  area. 

Dissimilar  exposure  and  soil  conditions  were  probably 
responsible  for  the  differences  between  the  windward 
communities  on  the  two  study  areas.  Limber  pine  was 
relatively  more  common  on  the  Zoratti  area  (t=4.391  p<0.001 
S  df=15;  Table  3)  and  a  greater  proportion  of  trees  of  both 
species  in  the  Zoratti  windward  community  exhibited  the 
prostrate  growth  form  characteristic  of  repeated  moisture 
stress  (Table  4:  Brown  1971c).  Exposure  explained  a 
significant  proportion  of  the  variation  in  the  DBH  (Table  5) 
of  Douglas  fir  and  the  height  (Table  6)  of  both  species  in 
the  Skelton  windward  community.  The  more  northerly  exposed 
plots,  although  classed  as  windward  on  the  basis  of  the 
associated  herbaceous  vegetation,  contained  trees  of  a  form 
more  typical  of  leeward  plots. 

Although  aspect  in  part  determines  the  intensity  of 
insolation  and  the  temperature  and  moisture  content  of  both 
the  soil  and  the  atmosphere  (Pearson  1920;  McMinn  1952),  it 
does  not  affect  plants  directly  (Billings  1952) .  Tree 
species  differ  in  their  ability  to  deal  with  the  available 
range  of  edaphic  and  atmospheric  conditions  so  that  exposure 
can  be  important  in  determining  the  distribution  and  local 
success  of  a  species.  Limber  pine  rely  upon  structural 
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Table  4. 


Proportion  of  prostrate  trees  in  the  windward 
communities  on  the  Zoratti  and  Skelton  areas. 


Study  Area 


Species 

Zoratti 

Skelton 

menziesii 

29.3% 

1.8% 

(27/92) 

(5/271) 

f lexilis 

30.5% 

4.2% 

(163/535) 

(5/117) 

.  ) 


( 
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Table  5.  Regression  statistics  describing  the  relationship 
between  diameter  at  breast  height  and  the  measured 
environmental  parameters.  Additional  statistics 
have  been  provided  in  Appendix  3. 


Independent  Variables 


Variablesc 1 > 

Beta <1  2 3  J 

R2 

not  in 

menziesii 

Windward 

Zoratti  {3} 
(25)  * 

Age  £4} 

”0.740 

0.548 

Slope 

Density 

Skelton  £2} 
(54) 

Leeward  Pure 

Age  £4} 
Exposure  £1} 
Density  £4} 

”0.583 

0.553 

0.276 

0.542 

Slope 

Zoratti  £2} 
(104) 

Age  £3) 
Density 

0.624 

”0.374 

0.633 

Slope 

Skelton  £2} 
(247) 

Leeward  Mixed 

Age  £3} 
Density  £4} 

0.699 

0.288 

0.570 

Slope 

Exposure 

Zoratti  £3} 
(28) 

Age  {2} 

0.823 

0.677 

Slope 

Density 

flexilis 

Windward 

Zoratti  {2} 
(94) 

Age  £2} 

Slope  £4} 

0.664 

0.210 

0.469 

Density 

Skelton  £4} 

(19) 

Leeward  Mixed 

Age  £4} 

0.652 

0.425 

Slope 

Density 

Exposure 

Zoratti 

(31) 

Age 

0.678 

0.459 

Slope 

Density 

1p<0.001  for  all  equations  except  Skelton  windward 
pine  0 .005>p>0 .00 1 . 

2Beta  (standardized  regression  coefficient)  is  a  measure  of 
the  average  effect  that  a  change  of  one  standard  deviation 
in  the  independent  variable  has  upon  the  dependent  variable 
(expressed  in  standard  deviations) . 

3Numbers  in  braces  refer  to  the  type  of  transformation 
applied:  {1}  squared;  £2}  square  root;  {3}  log;  {4}  inverse. 
Braces  following  the  study  area  name  indicate  that  the 
transformation  was  applied  to  the  dependent  variable. 

♦Sample  size. 
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Table  6.  Regression  statistics  describing  the  relationship 
between  height  and  the  measured  environmental 
parameters.  Additional  statistics  have  been 
provided  in  Appendix  3. 


Independent  Variables 

Variables*1*  Beta  Rz  not  in 


P .  menziesii 
Windward 


Zoratti  {3} 
(24) <3* 

*2*  Age  {4} 

-0.632 

0.399 

Density 

Slope 

Skelton  {2} 
(49) 

Exposure 
Age  {2} 

{1} 

0.467 

0.438 

0.304 

Slope 

Density 

Leeward  Pure 

Zoratti 

(88) 

Age  {3} 
Density 

0.667 

-0.277 

0.596 

Slope 

Skelton*  *  * 

Age  {3} 

0.701 

0.487 

Density 

(219) 

Exposure 

{4} 

-0.109 

Slope 

Leeward  Mixed 

Zoratti  {3} 
(25) 

Age  {3} 

0.792 

0.628 

Slope 

Density 

.  flexilis 

Windward 

Zoratti  {3} 
(83) 

Age  {3} 
Slope  {4} 

0.566 

0.255 

0.356 

Density 

Skelton  {3} 

(12) 

Exposure 

{1} 

0.720 

0.519 

Slope 

Density 

Age 

Leeward  Mixed 

Zoratti  {2} 
(26) 

Age  {2} 
Slope  {1} 

0.878 

-0.391 

0.551 

Density 

p<0.001  for  all 

equations  except 

Skelton 

windward 

pine  0.025>p>0.01 . 

2Numbers  in  braces  refer  to  the  type  of  transformation 
applied:  {1}  squared;  {2}  square  root;  {3}  log;  {4}  inverse. 
Braces  following  the  study  area  name  indicate  that  the 
transformation  was  applied  to  the  dependent  variable. 

3Sample  size. 

♦Height  of  Skelton  fir  pure  is  not  normally  distributed. 
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adaptations  such  as  stomatal  reduction,  thickened  epidermis, 
and  clustered  leaves  for  protection  against  water  loss,  but 
in  so  doing  restrict  their  ability  to  make  use  of  light  in 
photosynthesis.  Douglas  fir  and  spruce  have  adopted  an 
alternative  strategy,  reducing  water  loss  through 
physiological  adaptations  which  are  not  as  effective  in 
extreme  conditions  of  wind  and  dryness,  but  enhance 
photosynthesis  and,  therefore,  growth  (Bates  1923).  As  a 
result,  limber  pine  is  at  a  competitive  disadvantage  on  all 
but  the  more  xeric  sites  which  are  often  exposurally 
distinct  from  closed  canopy  Douglas  fir-spruce  forests 
(Pearson  1920;  Bates  1923;  Marr  1967). 

Although  increasing  slope  appears  to  have  had  a 
negative  effect  on  the  size  and  form  of  windward  limber  pine 
on  the  Zoratti  area  (Tables  5,  6,  and  7) ,  the  factor 
responsible  for  variation  was  probably  soil  depth  since  the 
slope  of  the  steepest  plots  was  dominated  by  the  inclination 
of  the  bedding  planes  of  the  exposed  bedrock.  Slope  may 
have  acted  more  directly  on  the  height  and  form  of  limber 
pine  in  the  leeward  mixed  forest.  This  community  appears  to 
be  advancing  slowly  on  to  flatter  ground  at  the  base  of  the 
ridge.  The  proportion  of  limber  pine  increases  with 
decreasing  slope  (Spearman* s  rho=-0.900  t=3.576  df=3 
0 . 05>p>0. 0 1)  suggesting  that  the  developing  pine  trees 
ameliorate  the  conditions  of  the  lower  slopes  allowing  the 
establishment  of  Douglas  fir  which  eventually  replace  the 
pine.  Reduced  light  conditions  in  the  presence  of  the 
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Table  7.  Regression  statistics  describing  the  relationship 
between  conformation  index  and  the  measured 
environmental  parameters.  Additional  statistics 


have  been  provided  in  Appendix  3. 


Independent  Variables 

Variables*1^  Beta  R2  not  in 


P.  menziesii 
Windward 
Zoratti 

(24) 

Skelton  {3} 
(49) 

Leeward  Pure 
Zoratti  {2} 
(38) 

Skelton  {3} 
(219) 

Leeward  Mixed 
Zoratti  {2} 

(25) 

P.  f lexilis 
Windward 
Zoratti 
(83) 

Skelton 

(12) 

Leeward  Mixed 
Zoratti 

(26) 


Age  {4} 

Slope  {4} 

-0.674 

0.298 

Density 
Age  {4} 

{1} 

-0.417 

-0.335 

Age  {4} 
Density 

{3} 

-0.354 

-0.312 

Age  {4} 
Density 

{3} 

-0.565 

-0.372 

Age 

0.829 

Age  {2} 

Slope  {1} 

0.670 

-0.200 

Density 

m 

-0.585 

Age 

Slope  {1} 
Density  {4} 

0.951 

-0.490 

0.338 

0.  573 

Density 

0.260 

Exposure 

Slope 

0.265 

Slope 

0.431 

Exposure 

Slope 

0.687 

Density 

Slope 

0.459 

Density 

0.343 

Age 

Exposure 

Slope 

0.647 


ip<0.001  for  all  equations  except  Skelton  windward 
fir  (0.005>p>0.001)  and  Skelton  windward  pine  (0. 05>p>0. 025) 
2Numbers  in  braces  refer  to  the  type  of  transformation 
applied:  {1}  squared;  {2}  square  root;  {3}  log;  {4}  inverse. 
Braces  following  the  study  area  name  indicate  that  the 
transformation  was  applied  to  the  dependent  variable. 

3Sample  size. 
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faster  growing  Douglas  fir  could  have  caused  limber  pine  on 
the  steeper  slopes  to  develop  longer,  relatively  thinner 
stems  (see  Brown  1971c). 

Variation  in  age  accounted  for  a  significant  proportion 
of  the  intraspecific  variation  in  the  DBH  (Table  5) ,  height 
(Table  6) ,  and  form  (Table  7)  of  conifers  in  the  three 
vegetation  types  examined,  with  the  exception  of  limber  pine 
in  the  Skelton  windward  community.  The  initial  attempt  at 
demonstrating  the  expected  correlation  between  the  age  and 
DBH  of  the  Skelton  windward  pine  was  unsuccessful  until 
three  of  the  small  sample  of  aged  trees  (21)  were  excluded 
from  the  analysis.  These  trees  were  probably  growing  at  a 
slower  than  normal  rate  since  two  forked  a  short  distance 
above  ground  level  and  the  third  had  partially  fallen  over. 
The  ineffectiveness  of  age  in  explaining  variation  in  the 
height  of  Skelton  pine  presumably  resulted  from  a  small 
sample  size  (n=12) ,  compounded  by  the  pattern  of  primary 
growth  of  limber  pine,  which  involves  a  loss  of  vigour  by 
the  original  leader  (Marr  1967)  and  production  of  a  more 
spherical  crown  than  is  typical  for  conifers  (see 
Assmann  1970,  pg.  39). 

Stand  density  was  another  important  factor  influencing 
the  size  (Tables  5  and  6)  and  form  (Table  7)  of  trees  in  the 
leeward  pure  forests  of  both  study  areas.  The  inverse 
relationship  between  density  and  DBH  and  height  is  a  well 
known  phenomenon  that  serves  as  the  basis  for  the 
silvicultural  practice  of  thinning  (Assmann  1970).  Crowding 
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restricts  the  size  of  the  crown  (Mitchell  1969)  and 
therefore  reduces  the  supply  of  photosynthetic  products 
required  for  the  growth  of  new  tissues  (Brown  1971a  and  b) . 

The  variation  in  age  and  density,  which  seemingly 
played  an  important  role  in  determining  the  form  of  the 
trees  available  to  porcupines,  was  probably  a  result  of  the 
pattern  of  forest  regeneration  following  burning.  Several 
limber  pine  and  Douglas  fir,  in  excess  of  200  years  of  age, 
were  present  on  both  study  areas.  The  densest  stands  in  the 
leeward  forests  generally  occurred  to  the  east  (downwind) 
of,  and  in  close  proximity  to,  these  older  trees,  suggesting 
that  stand  density  is  related  to  the  distance  from  a  seed 
source.  Trees  in  the  windward  community,  especially  limber 
pine,  appeared  to  be  relatively  more  successful  in  surviving 
fire,  presumably  having  been  protected  by  the  sparse  ground 
vegetation  and  the  open  nature  of  this  community.  Because 
of  a  more  evenly  distributed  seed  stock,  this  protection, 
however  slight,  may  explain  the  generally  greater  age  of  the 
windward  community  (Zoratti  t=2.228  0.05>p>0.02  S  df=243; 
Skelton  t=3.246  0.01>p  >0.001  S  df=109;  Table  3). 

Apparently  the  Zoratti  area  was  burned  somewhat  later  than 
the  Skelton  area  since  its  forests  were  appreciably  younger 
(t  =  4 . 2 09  p<0.001  S  df =550 ;  Table  3). 

At  the  time  of  measurement,  most  of  the  trees  on  the 
Skelton  area  were  older  and  larger  than  those  on  the  Zoratti 
area  (Tables  8,  9,  10a  and  11a).  Within  each  study  area, 
trees  in  the  leeward  forest  were  generally  larger  than 
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Table  8.  Descriptive  statistics  of  the  age  and  size  of 
individual  trees  including  the  mean,  standard 
deviation,  and  number  of  trees  examined. 


Age 

DBH 

Height 

(yr) 

(cm) 

(m) 

Cl 

P  •  menziesii 
Windward 


Zoratti  (X) 

59.5 

8.4 

3.7 

2.03 

(s) 

25.55 

6.18 

1.77 

0.837 

(n) 

26 

67 

65 

65 

Skelton 

74.2 

15.1 

6.6 

2.33 

18.36 

7.49 

2.93 

0.692 

54 

254 

235 

235 

Leeward  Pure 

Zoratti 

60.5 

17.1 

9.4 

1.67 

19.31 

9.92 

4.28 

0.592 

112 

378 

296 

296 

Skelton 

66.5 

19.9 

11.8 

1.61 

16.08 

10.47 

4.25 

0.541 

246 

1218 

1057 

1057 

Leeward  Mixed 

Zoratti 

47.4 

11.4 

5.4 

1.  74 

19.83 

9.25 

3.46 

0.715 

27 

67 

62 

62 

P.  flexilis 

Windward 

Zoratti 

68.  8 

9.  1 

3.  1 

2.70 

26.78 

6.12 

1.27 

1.042 

1  1 1 

370 

299 

299 

Skelton 

74.9 

12.6 

4.3 

2.91 

21.93 

7.35 

1.75 

0.  930 

22 

112 

96 

96 

Leeward  Mixed 

Zoratti 

50.6 

11.5 

4.8 

2.  12 

17.  92 

6.  76 

1.96 

0.  831 

28 

94 

75 

75 

Picea  spp. 

Leeward  Pure 

Zoratti 

63.2 

24.7 

12.4 

2.00 

15.47 

11.26 

4.3 

0.  641 

8 

49 

46 

46 
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Table  9.  Comparisons  of  the  ages  of  conifers  on  the  Zoratti 
and  Skelton  areas.  Comparison  of  Douglas  fir  in 
the  three  Zoratti  communities  is  presented  in 
Table  13. 

Mean 


Comparison1 

(yr) 

t 

P 

df 

Type 

ZWF 

59.5 

1.61 

p>0. 1 

135 

P 

ZPW 

68.8 

SWF 

74.2 

0.14 

p>0. 75 

74 

P 

SWP 

74.9 

ZWF 

59.5 

2.63 

0 . 02>p>0 .01 

38 

S 

SWF 

74.2 

ZWP 

68.8 

1.00 

p>0 .25 

131 

P 

SWP 

74.9 

ZPF 

60.5 

2.87 

0.01>p>0. 001 

184 

S 

SPF 

66.5 

ZPF 

60.5 

0.39 

p>0. 5 

118 

P 

ZPS 

63.2 

ZMF 

47.4 

0.64 

p>0.5 

53 

P 

ZMP 

50.6 

SWF 

74.2 

3.08 

0.01>p>0. 001 

298 

P 

SPF 

66.5 

ZWP 

68.8 

4.28 

p<0.001 

61 

S 

ZMP 

50.6 

omparison  code: 

First  character 

z-: 

Zoratti 

S-Skelton 

Second  character  W-i 

Windward 

P-Leeward  pure  M 

-Leeward 

Third  character 

F-j 

Douglas 

fir  P-Limber 

pine 

S-S  pruce 

c 
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Table  10.  Comparisons  of  the  diameter  at  breast  height  of 
trees  on  the  Zoratti  and  Skelton  study  areas. 
Comparison  of  Douglas  fir  in  the  three  Zoratti 
communities  is  presented  in  Table  13.  The 
extreme  right-hand  column  bears  the  following  code 
for  interpretation  of  the  analyses  of  covariance: 

1  -  Similar  growth  rates  throughout; 

2  -  Different  initial  growth  rates; 

3  -  Different  growth  rates  throughout, 

a)  Difference  in  growth  curves  corrected  by 
log-log  transformation 

b)  Interaction  -  the  slopes  of  the  individual 


growth  curves  are  not  equivalent. 


a)t-test  Adjusted  b)  Analysis  of  Covariance 
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3This  F-test  examines  the  contribution  of  age  in  explaining  variation  in  DBH. 
♦Explanation  of  this  comparison  code  is  provided  in  Table  9. 

5The  dependent  variable  has  been  log  transformed  to  correct  for  mildly  heterogeneous 
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Table  11.  Comparisons  of  the  height  of  trees  on  the  Zoratti 
and  Skelton  study  areas.  Comparison  of  Douglas 
fir  in  the  three  Zoratti  communities  is 
presented  in  Table  13.  The  extreme  right-hand 
column  bears  the  following  code  for  interpretation 
of  the  analyses  of  covariance: 

1  -  Similar  growth  rates  throughout; 

2  -  Different  initial  growth  rates; 

3  -  Different  growth  rates  throughout, 

a)  Difference  in  growth  curves  corrected  by 
log-log  transformation 

b)  Interaction  -  the  slopes  of  the  individual 
growth  curves  are  not  equivalent. 


a)t-test  Adjusted  b)  Analysis  of  Covariance 
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windward  trees,  even  though  the  windward  forests  were  older. 
Although  fir  tended  to  be  larger  than  pine  within  the 
windward  community,  there  was  no  significant  difference  in 
the  average  size  of  these  two  species  in  the  leeward  mixed 
forest.  Finally,  the  spruce  in  the  leeward  pure  forest  on 
the  Zoratti  area  were  larger,  though  not  older,  than  the 
surrounding  Douglas  fir. 

Conformation  of  trees  also  varied  between  forest  types 
and  study  areas  (Tables  12a  and  13).  (Interspecific 
comparisons  involving  limber  pine  do  not  appear  in  Table  12 
because  of  its  unusual  pattern  of  primary  growth.)  The 
crowns  of  windward  trees  and  leeward  trees  in  low  density 
stands  often  extended  to  within  a  meter  of  the  ground, 
whereas  often  only  the  upper  third  of  leeward  trees  in  dense 
stands  was  foliated.  Windward  trees  had  stout  trunks, 
compared  to  leeward  trees  of  similar  height,  probably 
because  of  the  stimulatory  effect  of  wind-induced  sway  on 
secondary  growth  (Jacobs  1954) .  The  Douglas  fir  in  the 
leeward  pure  forest  on  the  Zoratti  area  had  more  extensive 
crowns  on  average  than  their  counterparts  on  the  Skelton 
area  as  a  result  of  a  higher  mean  density  on  the  latter  area 
(Zoratti  weighted  mean  (±  se)  12.18  ±  0.286  trees/100  sg.  m 
(n=359)  ;  Skelton  weighted  mean  (±  se) 

17.08  ±  0.163  trees/100  sg.  m  (n=1221);  t=5.61  p<0.001  P 
df =1578) . 

The  observed  similarities  and  differences  in  the  size 
and  form  of  trees  can  be  interpreted  by  three  complementary 
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Table  12.  Comparisons  of  the  conformation  index  of  trees 
on  the  Zoratti  and  Skelton  study  areas. 

Comparison  of  Douglas  fir  in  the  three  Zoratti 
communities  is  presented  in  Table  13.  The 
extreme  right-hand  column  bears  the  following  code 
for  interpretation  of  the  analyses  of  covariance: 

1  -  Similar  growth  rates  throughout; 

2  -  Different  initial  growth  rates; 

3  -  Different  growth  rates  throughout, 

a)  Difference  in  growth  curves  corrected  by 
log-log  transformation 

b)  Interaction  -  the  slopes  of  the  individual 
growth  curves  are  not  equivalent. 


a)t-test  Adjusted  b)  Analysis  of  Covariance 
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Table  13.  Comparisons  of  the  age,  diameter  at  breast  height, 
height  and  conformation  index  of  Douglas  fir  in 
the  three  conifer  dominated  communities  on  the 
Zoratti  study  area. 

Age  -  One-way  analysis  of  variance 

F=4 . 523  0. 01>p>0 .025  2  and  162  df 

Mixed  Windward  Pure 

47.4  59.  5  _ 60. 5  (Student-Newman-Keuls  test  p=0.05) 

log  DBH 

a)  One-way  analysis  of  variance 

F=3 1.790  p<0.001  2  and  509  df 

Windward  Mixed  Pure 

8.4 _ 11.4  17.1  (Student-Newman-Keuls  test  p=0.05) 

b)  Analysis  of  covariance 

covariate  -  log  Age  F=177.281  p<0.001 

F=1 1 . 457  p<0.001  2  and  161  df 

Windward  Pure  Mixed  (adjusted  means) 

6.6  12.3  13.2  (Student-Newman-Keuls  test  p=0.05) 

log  Height 

a)  One-way  analysis  of  variance 

F=75. 822  p<0.001  2  and  420  df 

Windward  Mixed  Pure 

3.7  5.4  9.4  (Student-Newman-Keuls  test  p=0.05) 

b)  Analysis  of  covariance 
significant  interaction 

F=9.607  p<0.001  2  and  136  df 

Cl  -  Heterogeneous  variances 

a)  Approximate  F-test  (Sokal  and  Rohlf  1969:  376-379) 

Fs=5. 238  0. 01 >p>0. 005  2  and  102  df 

Pure  Mixed  Windward 

1^6 7 _ 1.74  2.03  (based  on  t-tests  p=0.05) 
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growth  rate  relationships  revealed  by  analyses  of  covariance 
(Tables  10b,  11br  12b  and  13),  In  applying  an  analysis  of 
covariance,  several  assumptions  must  be  satisfied.  If  the 
variances  under  examination  are  not  homogeneous,  even 
following  a  suitable  transformation,  this  test  is 
inappropriate  (Steel  and  Torrie  1960).  Additionally,  the 
linear  regression  equations  explaining  the  association 
between  the  independent  variable  (age)  and  the  dependent 
variable  (DBH,  height,  or  Cl)  for  each  of  the  samples  must 
share  a  common  slope.  Failure  to  fulfill  this  requirement 
also  invalidates  the  analysis  of  covariance,  but  in  the 
present  context  the  demonstrated  interaction  (i.e.  different 
slopes)  signifies  different  growth  rates  for  the  two  groups 
of  trees  being  compared. 

The  inability  to  demonstrate  a  difference  between  two 
means  after  adjusting  for  the  affect  of  age  suggests  similar 
growth  rates  during  all  periods  of  development.  If  these 
means  differed  prior  to  adjustment  for  the  effect  of  the 
independent  variable  (as  demonstrated  by  a  t-test) ,  the 
inequality  was  probably  due  to  a  difference  in  age.  A 
significant  difference  between  the  adjusted  means  may  arise 
for  two  different  reasons.  If  the  dependent  variable 
remains  untransformed,  disparity  between  means  reveals 
dissimilar  Y-intercepts  for  the  regression  equations  (see 
Figure  5a  and  b) •  Since  tree  growth  tends  to  follow  a 
sigmoid  curve  (Husch  et  al.  1972)  and  the  majority  of  trees 
were  measured  after  they  were  greater  than  fifty  years  of 
age,  extrapolation  of  the  growth  curves  to  age  zero  is 
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Figure  5.  The  effect  of  log  transformation  on  growth  curves. 

A  and  B  are  drawn  from  the  regression  equation 
used  in  the  analysis  of  covariance  comparing  the 
height  of  Douglas  fir  trees  in  the  leeward  pure 
forests  on  the  two  study  areas  (Zoratti  a  : 

Skelton  ©  ) .  C  and  D  are  drawn  from  similar 
equations  used  in  comparing  the  height  of  Douglas 
fir  in  the  windward  communities  iZoratti  ❖  : 
Skelton  a  ) • 
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unrealistic.  The  observed  difference  therefore  indicates 
unequal  initial  growth  rates.  Conversely,  if  both  the 
dependent  and  independent  variables  were  log  transformed 
prior  to  analysis,  the  trees  in  question  must  have  grown  at 
perennially  different  rates  (Figure  5c  and  d)  . 

Douglas  fir  in  both  the  windward  and  leeward 
communities  on  the  Skelton  area  enjoyed  superior  rates  of 
height  growth  (at  least  initially)  relative  to  fir  on  the 
Zoratti  area,  even  though  there  was  no  difference  between 
areas  in  the  rates  of  diameter  growth.  While  Douglas  fir 
grew  more  rapidly  in  diameter  than  pine  in  the  windward 
communities,  the  interspecific  difference  in  the  leeward 
mixed  community  was  not  statistically  significant.  On  both 
study  areas,  growth  appears  to  have  been  most  vigorous  on 
leeward  sites.  The  comparison  of  spruce  and  Douglas  fir  in 
the  Zoratti  leeward  pure  forest  is  unique  in  that  spruce 
grew  significantly  faster  in  both  DBH  and  height  and  also 
maintained  a  more  rapid  increase  in  conformation  index. 

Most  of  the  transformations  required  by  the  analyses  of 
covariance  suggest  that  the  trees  on  both  study  areas,  with 
the  possible  exception  of  those  in  the  Zoratti  leeward  mixed 
sites,  are  beyond  the  period  of  most  rapid  growth.  Even  on 
poor  sites,  interior  Douglas  fir  reaches  the  point  of 
inflection  of  the  growth  curve  at  about  40  years  of  age 
(Forestry  Handbook  for  British  Columbia  1971),  after  which 
the  size  of  annual  increments  begins  to  decrease  (Husch 
et  al.  1972) . 


Sq 


' 


. 

- 


50 


Characteristics  of  Trees  Used  by,  Porcupines 

The  relative  proportions  of  used  trees,  as  revealed  by 
the  presence  of  feeding  scars,  are  an  expression  of  the 
interspecific  and  interhabitat  preferences  of  porcupines. 
Analysis  of  this  information  was  based  on  a  sample 
representing  many  years  of  feeding,  so  that  the  influence  of 
individual  preferences  is  suppressed  and  the  resulting 
conclusions  reflect  only  the  majority  response.  These 
preferences  are  assumed  to  be  constant  through  time. 

Interspecific  Preferences 

Interspecific  comparisons  of  the  proportion  of  used 
trees  suggest  that  porcupines  that  have  inhabited  the  two 
study  areas  did  not  prefer  one  conifer  species  over  another 
within  a  given  community  (Tables  14  and  15).  This 
observation  conflicts  with  the  published  accounts  of 
porcupine  food  habits,  which  generally  characterize  this 
species  as  a  selective  feeder.  The  accuracy  of  many  of 
these  reports  is,  however,  difficult  to  assess  because  data 
were  not  provided  regarding  the  availability  of  the  tree 
species  involved  and  preference  rankings  are  apparently 
subjective  (Gabrielson  1928;  Taylor  1935;  Curtis  and 
Kozicky  1944;  Shapiro  1949).  Two  other  studies  describe 
interspecific  preferences  that  are  unsupported  by  the  data 
presented.  Gill  and  Cordes  (1972)  stated  that  limber  pine 
was  the  preferred  species  in  low  elevation  krummholz, 
similar  to  the  windward  communities  described  here,  even 
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Table  14.  Proportion  of  trees  used  within  species  and 

community  categories  from  the  Zoratti  and  Skelton 
study  areas. 


Number 

Number  Not  % 

Used  Used  Used 


Windward 


Zoratti 


P.  menziesii 

11 

56 

16.4 

P.  flexilis 

96 

274 

25.9 

Total 

107 

330 

24.5 

Skelton 

P.  menziesii 

37 

217 

14.6 

P.  flexilis 

9 

104 

8.0 

Total 

46 

321 

12.5 

Leeward  Pure 

Zoratti 

P.  menziesii 
High  Density 

119 

76 

61.0 

Low  density 

124 

59 

67.8 

Total 

243 

135 

64.3 

Picea  spp.  1 

36 

13 

73.5 

Skelton 

P.  menziesii 
High  Density 

147 

7  02 

17.3 

Low  Density 

164 

205 

44.4 

Total 

31  1 

9C7 

25.5 

Leeward  Mixed 
Zoratti 


P.  menziesii 

23 

44 

34.3 

P.  flexilis 

22 

72 

23.  4 

Total 

45 

116 

28.0 

4Spruce  and  Douglas  fir  were  not  sampled  with  the 
same  techniques  (see  Methods) ,  consequently  spruce 
does  not  participate  in  the  Zoratti  leeward  pure 
total. 
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Table  15,  Selected  inter-  and  intraspecific  comparisons  of 
the  proportion  of  trees  used  by  porcupines  on 
the  two  study  areas.  Each  entry  includes  the 
log  likelihood  ratio  statistic  (G)  and  the 
probability  of  a  similar  or  greater  difference 
between  groups.  One  degree  of  freedom  is  associated 
with  each  test. 


ZWP* 


r~ — — — ■ — 

G=2. 448 

p>0  •  1 


ZPF 


ZPS 


ZMF 


ZMP 


SWF 


G=53.203 
p<0. 005 


G=4 • 836 
0 . 05>p 
p>0 .025 


G=0. 034 
p>0 .5 


— ^ - — 

G=1. 276 


p>0. 1 


G= 1 9 . 649 
p<0 . 005 


■ - - - T~~ - — - 1 

G=0. 1 40  |  G=1 .793 
p>0. 5  |  p>0. 1 


— — — f 

G= 17.970| 

I 

p<0.005  | 

- < 


SWP 


SPF 


G=182.568 
p<0 . 005 


ZWF 


ZPF 


ZWP 


ZMF 


G=2. 706  j 

i 

P=0.1  I 
~ - - — H 

G= 14. 584 | 

l 

p<0.005  | 


SWF 


lConsult  Table  9  for  an  explanation  of  the  comparison  codes. 
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though  it  was  the  only  species  available.  Krefting  et  al. 
(1962)  ranked  the  four  major  tree  species  on  their  Wisconsin 
study  area  according  to  porcupine  use  but  re-analysis  of 
these  data  indicate  that  the  proportion  of  used  trees  does 
not  vary  statistically  between  species  (G=1.634  df=3  p>0.5). 
Although  porcupines  in  some  areas  do  exhibit  interspecific 
preferences  (Curtis  1941  ;  Rudolf  1949;  Brander  1973)  this 
habit  is  apparently  not  universal. 

Individual  trees  are  occasionally  used  during  several 
years  (see  Tables  25  and  26) .  If  these  trees  are  weighted 
for  the  number  of  times  they  have  been  used,  i.e.  a  tree 
used  in  two  different  years  is  recorded  as  two  used  trees, 
spruce  were  used  proportionately  more  often  than  fir  in  the 
Zoratti  leeward  pure  forest  (G=9.943  df = 1  p<0.005)  and  fir 
were  used  more  frequently  than  pine  in  the  Skelton  windward 
community  (G=5.367  df=1  0 . 025>p>0.01 ) .  However,  it  was  not 
possible  to  weight  trees  for  the  number  of  times  they  were 
rejected  by  porcupines  so  that  unused  trees  may  be 
underrepresented. 

Intercommunity  Preferences 

In  contrast  to  this  apparent  lack  of  interspecific 
discrimination,  porcupines  on  the  Skelton  and  Zoratti  areas 
exhibited  marked  community  preferences  (Tables  14,  15, 
and  16) .  Within  a  study  area,  the  communities  most  amenable 
to  tree  growth  received  the  greatest  attention  from 
porcupines.  Douglas  fir,  in  both  the  leeward  pure  and 
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Table  16.  Selected  intercommunity  comparisons  of  the 

proportion  of  trees  used  by  porcupines  on  the  two 
study  areas.  Each  entry  includes  the  log 
likelihood  ratio  statistic  (G)  and  the  probability 
of  a  similar  or  greater  difference  between  groups. 
One  degree  of  freedom  is  associated  with  each 
test. 


Zoratti 

Leeward 

Pure 

Zoratti 

Leeward 

Mixed 

Skelton 
Win dward 


Skelton 

Leeward 


rG=132. 637 

p<0 . 005 


G=0 . 567 


p>0  •  1 


G=  18. 224 


p<0. 005 


G=59. 639 


p<0. 005 


-f 


G=1 82 • 56  8 
p<0. 005 


G=29.210 | 

I 

p<0.005  | 

— —j 


Zoratti  Zoratti  Skelton 

Windward  Leeward  Windward 

Pure 
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leeward  mixed  forests,  were  used  proportionately  more 
frequently  than  windward  fir.  The  seeming  preference  for 
leeward  pure  fir  over  fir  in  the  leeward  mixed  forests  on 
the  Zoratti  area  probably  reflects  the  relative  lengths  of 
time  these  two  communities  have  been  subject  to  porcupine 
feeding  rather  than  preferential  utilization.  As  a  whole, 
the  Zoratti  leeward  mixed  community  was  used  in 
approximately  equal  proportion  to  the  windward  community 
because  of  an  equivalent  level  of  feeding  on  limber  pine. 

The  relationship  between  use  by  porcupines  and 
individual  tree  vigour  has  been  widely  documented 
(Rudolf  1949;  Curtis  and  Wilson  1953;  Krefting  et  al.  1962) 
and  will  be  examined  in  more  detail  later.  If  the 
porcupines  that  have  fed  on  the  Zoratti  and  Skelton  areas 
selected  trees  on  the  basis  of  vigour  the  collective 
response  would  be  manifest  in  intercommunity  preferences 
similar  to  those  observed.  Lack  of  interspecific 
preferences,  in  the  light  of  differences  in  growth  rates, 
suggests  that  the  importance  of  vigour  is  not  consistent 
between  species.  The  influence  of  the  intensive  winds 
common  in  the  windward  communities  on  the  selection  of 
feeding  trees  by  porcupines  is  uncertain,  since  Douglas  fir 
in  this  community  on  the  Zoratti  area  were  used  less 
frequently  than  leeward  mixed  fir,  but  limber  pine  in  these 
two  forest  types  were  used  in  equivalent  proportions. 

Affirmation  of  intercommunity  preferences  is  provided 
by  the  outcome  of  multiple  regression  analyses  relating  the 
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proportion  of  trees  used  by  porcupines  per  plot  to  the 
proportion  of  Douglas  fir  per  plot,  plot  density,  slope, 
exposure,  mean  DBH,  height,  conformation  index  and  age.  Of 
these,  the  two  variables  that  were  the  most  effective 
predictors  of  the  intensity  of  use  on  the  two  study  areas 
(Table  17)  also  allowed  the  greatest  discrimination  between 
communities  (Zoratti  -  proportion  of  Douglas  fir/plot; 
Skelton  -  exposure) .  Variables  that  subsequently  entered 
the  equations  were  more  useful  in  differentiating  between 
stands  within  a  vegetation  type  than  between  types. 

Stand  density  is  the  most  important  single  variable  in 
explaining  variation  in  the  proportion  of  trees  used  within 
a  plot,  but  because  density  is  influential  in  determining 
the  size,  form  and  growth  rate  on  both  study  areas,  its  role 
in  food  selection  by  porcupines  may  be  multifaceted. 

Density  may  have  been  considered  directly  by  porcupines 
feeding  in  at  least  the  Zoratti  windward  community,  as  the 
dispersion  of  trees  did  not  account  for  a  significant 
proportion  of  variation  in  the  size  or  form  of  either  tree 
species  (Tables  5,  6  and  7) .  The  inverse  relationship 
between  density  and  the  intensity  of  use,  also  reported  by 
van  Deusen  and  Myers  (1962) ,  may  not  be  consistent 
throughout  the  range  of  possible  densities,  since  Curtis  and 
Wilson  (1953)  found  that  use  increased  with  density  in  very 
dispersed  stands  of  Ponderosa  pine  (Pinus  ponderosa:  mean 
density=3.32  trees/100  sg.  m) . 


, 

> 


.  ( 

■  V 

' 

. 

( 

.  {'  ■  , 


\  ...  ! .  :  fcjii  frill  d  ••  )  :\i . 


57 


Table  17.  Regression  statistics  describing  the  relationship 
between  the  proportion  of  trees  used  by  porcupines 
per  plot  and  various  plot  attributes.  Additional 
statistics  have  been  provided  in  Appendix  3. 


Independent 

Variables  Beta  R1 2  p 


Windward 


Zoratti 
(9)  <2-» 

Density  {4} <  *  * 

0.  744 

0.  553 

0. 025>p>0 .01 

Skelton 

HI) 

p>0 . 05 

Leeward  Pure 
Zoratti  {5} 
(12) 

Density  {1}  <  3  * 
Slope  £4} 

-0.833 

-0.690 

0.499 

0 . 05>p>0 .025 

Skelton 

(29) 

Density  £3} 
Slope  £3} 
Exposure  £1} 

-0.500 

0.396 

0.307 

0.717 

p<0.001 

Leeward  Mixed 
Zoratti 
(5) 

#f ir/#total 

0.958 

0.917 

0. 025>p>0. 01 

Combined 

Zoratti 

(27) 

#f ir/#total 
Mean  DBH  £2} 

0.456 

0.429 

0.713 

p<0 . 001 

Skelton 

t40) 

Exposure  £1}<4*0.479 
Density  £4}  0.492 

0.628 

pCO.001 

Slope  {4}  -0.235 


1Numbers  in  braces  refer  to  the  type  of  transformation 
applied:  {1}  squared;  {2}  square  root;  {3}  log;  {4}  inverse, 
{5}  arcsin.  Braces  following  the  study  area  name  indicate 
that  the  transformation  was  applied  to  the  dependent 
variable. 

2Sample  size. 

Proportion  of  Douglas  fir  was  the  first  independent 
variable  to  enter  this  equation,  but  it  dropped  out  after 
the  addition  of  slope. 

4Mean  DBH  was  the  first  independent  variable  to  enter 
this  equation,  but  it  dropped  out  after  the  addition  of 
density. 
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The  role  of  slope  in  the  selection  of  stands  in  the 
leeward  pure  forests  is  unclear  because  slope  is  apparently 
not  important  in  explaining  variation  in  tree  stature 
(Tables  5,  6  and  7).  Slope  could  be  related  to  an  aspect  of 
porcupine  biology  other  than  feeding,  or  it  could  be  highly 
correlated  with  unmeasured  factors  that  influence  porcupines 
more  directly.  Although  van  Deusen  and  Myers  (1962) 
reported  a  similar  correlation  between  slope  and  the 
proportion  of  used  trees,  they  likewise  offered  no 
explanation  of  its  importance. 

The  regression  analysis  for  the  Zoratti  leeward  mixed 
forest  is  presented  with  reservation.  Although  the 
proportion  of  Douglas  fir  per  plot  was  the  only  variable  to 
explain  a  significant  portion  of  the  variation  in  the 
intensity  of  use,  a  thorough  analysis  was  not  possible 
because  the  number  of  variables  in  the  equation  rapidly 
approached  the  number  of  observations  (Draper  and 
Smith  1966) .  The  proportion  of  fir  per  plot  was  also  the 
first  variable  to  enter  the  regression  equation  in  the 
Zoratti  leeward  pure  analysis,  but  it  was  removed  from  the 
equation  following  the  addition  of  density  and  slope  because 
it  had  ceased  to  make  a  significant  contribution  to  the 
explanation  of  variation. 

Study  Area  Differences 

The  Zoratti  area  has  been  more  extensively  used  by 
porcupines  (Tables  14  and  16)  than  the  Skelton  area,  but 
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because  individual  porcupines  probably  could  not  have 
exercised  a  choice  between  these  two  areas,  inferences 
regarding  the  relative  quality  of  these  forests  would 
probably  be  unfounded.  Taylor  *1935)  stated  that  the 
prevalence  of  porcupine  feeding  in  a  particular  forest  is 
primarily  determined  by  the  number  of  porcupines  present, 
which  is  greatly  affected  by  the  physical  location  of  the 
forest  and  need  not  be  related  tc  the  type  of  trees 
available.  Both  Taylor  (1935)  and  van  Deusen  and 
Myers  (1962)  have  reported  that  trees  on  the  periphery  of  a 
forest  receive  the  greatest  attention  from  porcupines  and 
that  this  habit  was  related  to  the  proximity  of  the  forest 
edge  to  summer  feeding  areas.  The  Zoratti  area  is  a  narrow, 
linear  ridge  parallelling  the  edge  of  the  prairie  so  that 
the  entire  forest  is  essentially  peripheral,  whereas  the 
Skelton  area  lies  perpendicular  to  the  prairie  and  is  a 
larger,  more  continuous  block  of  timber.  This  difference  in 
size  and  orientation  may  explain  the  difference  in 
utilization  between  the  two  areas. 

Individual  Tree  Preferences  -  Relative 

The  preferential  use  of  a  specific  community  or  type  of 
stand  by  porcupines  ultimately  depends  upon  the  deliberate 
use  of  individual  trees.  Comparison  of  the  size  and 
conformation  of  used  trees  relative  to  unused  trees  could 
disclose  the  existence  and  form  of  porcupine  preferences. 

The  average  position  of  each  tree  in  relation  to  the 
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community  mean,  during  the  20  years  prior  to  measurement, 
has  been  examined  rather  than  the  tree9s  absolute  size  or 
age  because,  with  the  exception  of  the  leeward  pure  forests, 
only  a  small  proportion  of  the  available  trees  was  used  by 
porcupines  during  a  given  year.  Z-transf ormations  (Xi-X/s; 
Sokal  and  Rohlf  1969)  of  age,  DBH,  height  and  conformation 
index  were  calculated  for  each  tree  for  each  of  the  20  years 
prior  to  measurement.  The  average  transformation  of  the 
four  variables  for  each  tree  during  this  period  is  a  measure 
of  its  relative  status  in  the  community,  measured  in 
standard  deviations  from  the  mean  of  zero. 

Porcupines  have  used  trees  that  were  commonly  larger 
and  had  more  extensive  crowns  than  the  trees  they  neglected 
(Tables  18a,  19a  and  20a) ,  even  though  these  trees  were 
generally  not  older  than  unused  trees  (Table  21).  Although 
analyses  of  covariance  could  not  be  applied  universally,  the 
results  indicate  that  the  trees  used  by  porcupines  were 
faster  growing  than  unused  trees  (Tables  18b,  19b  and  20b) , 
an  observation  which  has  been  widely  reported  (Rudolf  1949; 
Curtis  and  Wilson  1953;  Krefting  et  al.  1962;  Spencer  1964). 
Discriminant  analyses  (Table  22)  indicate  that  DBH  and 
growth  form  are  more  effective  than  age  or  height  at 
distinguishing  used  trees  from  unused  trees  and  suggest  that 
any  tendency  for  used  trees  to  be  older  or  taller  was  an 
artifact  of  the  strong  intercorrelations  between  all  four 
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Table  18,  Comparisons  of  the  relative  diameters  at  breast 
height  of  porcupine  feeding  trees  with  unused 
trees  on  the  Zoratti  and  Skelton  areas.  The 
sizes  of  used  and  unused  limber  pine  in  the 
Skelton  windward  community  were  not  compared 
because  only  two  trees  were  used  during  the 
20  years  preceding  measurement. 


a)t-test  Adjusted  b)  Analysis  of  Covariance 
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Table  19.  Comparisons  of  the  relative  heights  of  porcupine 
feeding  trees  with  unused  trees  on  the  Zoratti 
and  Skelton  areas.  The  sizes  of  used  and  unused 
limber  pine  in  the  Skelton  windward  commumity 
were  not  compared  because  only  two  trees  were 
used  during  the  20  years  preceding  measurement. 


a)t-test  Adjusted  b)  Analysis  of  Covariance 
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Table  20.  Comparisons  of  the  relative  conformation 

indices  of  porcupine  feeding  trees  with  unused 
trees  on  the  Zoratti  and  Skelton  areas.  The 
form  of  used  and  unused  limber  pine  in  the 
Skelton  windward  community  was  not  compared 
because  only  two  trees  were  used  during  the 
20  years  preceding  measurement. 


a) t-test  Adjusted  b)  Analysis  of  Covariance 
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Table  21.  Comparisons  of  the  relative  ages  of  porcupine 
feeding  trees  with  unused  trees  on  the  Zoratti 
and  Skelton  areas.  The  ages  of  used  and  unused 
limber  pine  in  the  Skelton  windward  commumity 
were  not  compared  because  only  two  trees  were 
used  during  the  20  years  preceding  measurement. 


Means 

t 

P 

df 

Type 

ZWF1 used 

0.26 

0.  81 

p>0  •  4 

23 

P 

unused 

-0.  10 

ZPF  used 

0.05 

0.58 

p>0. 5 

87 

S 

unused 

-0.06 

ZMF  used 

0.  73 

2.02 

p>0 .05 

26 

P 

unused 

-0.24 

ZWP  used 

-0.26 

1.61 

p>0 . 1 

38 

S 

unused 

0.05 

ZMP  used 

1.05 

2.  46 

0. 05>p>0 . 02 

28 

P 

unused 

-0.19 

ZPS  used 

-0.17 

0.  80 

p>0  •  4 

6 

P 

unused 

0.50 

SWF  used 

0.19 

0,  44 

p>0.5 

51 

P 

unused 

-0.02 

SPF  used 

0.40 

2.98 

0. 01>p>0.001 

43 

s 

unused 

-0.06 

Explanation  of  this  comparison  code  is  provided 
in  Table  9. 
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Table  22.  Discriminant  analysis  statistics  detailing  the 

measured  variables  that  best  distinguished  between 
porcupine  feeding  trees  and  unused  trees  on  the 
Zoratti  and  Skelton  areas.1 


Independent  Discriminant  Wilks*  Chi-square 

Variable2  Constant  Coefficient  Lambda  (1  df) 
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In  addition  to  increased  shoot  and  xylem  growth, 
vigorous  trees  produce  a  larger  annual  increment  of  phloem 
than  suppressed  trees  (Bannan  1955;  Grillos  and  Smith  1959) . 
Rapidly  growing  trees  with  large  open  crowns  therefore 
provide  more  potential  food  for  porcupines  because  they 
support  more  foliage  and  a  greater  volume  of  inner  bark  than 
crowded  trees.  Porcupines  generally  feed  on  inner  bark  that 
is  accessible  from  positions  that  can  be  maintained  without 
undue  exertion  (Taylor  1935;  Spencer  1964).  This  type  of 
tree,  which  is  characterized  by  more  large  branches,  would 
therefore  allow  porcupines  greater  access  to  the  potential 
food. 

Diameter  at  breast  height  and  conformation  index  both 
incorporate  information  about  crown  width  and  stem  diameter 
and  relate  a  measure  of  the  potential  food  in  a  tree  and  its 
accessibility,  A  porcupine  that  was  cognizant  of  this 
relationship  could  therefore  assess  at  least  the  quantity  of 
food  available  in  a  given  tree  without  having  to  expend 
energy  climbing  it.  Presumably,  trees  can  become  too  large 
to  be  successfully  climbed  (Taylor  1935;  Curtis  and 
Kozicky  1944;  Curtis  and  Wilson  1953),  resulting  in  a 
decline  in  their  use  by  porcupines  even  though  they  contain 
large  quantities  cf  high  quality  food. 

Individual  Tree  Preferences  -  Absolute 

Only  trees  in  the  leeward  pure  forests  were  fed  upon 
with  sufficient  frequency  to  allow  an  examination  of  the 


' 

. 

1 1 


- 


71 


absolute  size  of  used  trees.  The  tendency  for  porcupines  to 
use  trees  that  were  larger  than  unused  trees  (Figures  6-11) 
has  already  been  noted.  Although  the  mean  DBH,  height  and 
growth  form  of  used  trees  in  the  Zoratti  forest  increased 
through  the  20  years  prior  to  measurement,  the  average  DBH 
and  height  in  the  Skelton  forest  remained  relatively 
constant  (Table  23  and  Figures  6-11). 

The  average  diameters  of  used  trees  from  both  study 
areas  fall  within  the  range  of  size  preferences  reported 
from  various  regions  of  North  America  (Taylor  1935;  Curtis 
and  Wilson  1953;  Krefting  et  al.  1962;  van  Deusen  and 
Myers  1962) .  The  predilection  of  porcupines  for  a  specific 
size  class  of  tree  irrespective  of  geographic  location  or 
the  species  of  tree  involved,  lends  support  to  the 
suggestion  of  Curtis  and  Wilson  (1953)  that  porcupines  may 
be  most  adept  at  climbing  trees  between  15  and  25  cm  DBH;  a 
physical  limitation  which  could  have  several  consequences  on 
the  pattern  of  forest  use.  If  porcupines  favoured  a 
particular  size  of  tree,  the  scope  for  additional 
discrimination  on  the  basis  of  other  attributes,  such  as 
vigour,  would  be  reduced.  In  a  relatively  even-aged  forest, 
vigorous  trees  of  the  appropriate  dimensions  would  only  be 
available  for  a  short  period  of  time  relative  to  the 
existence  of  the  forest  and  the  intensity  of  porcupine  use 
would  be  expected  to  peak  and  then  decline  as  vigorous  trees 
became  too  large  to  climb. 
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Figure  6.  Means  and  standard  errors  of  the  diameters  at 

breast  height  of  Douglas  fir  used  (right  bar 
of  each  pair)  and  not  used  by  porcupines  in 
the  leeward  pure  forest  on  the  Zoratti  study 
area  for  each  of  the  20  years  prior  to 
measurement.  Also  included  for  each  year  are 
the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t-test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Figure  7. 


Means  and  standard  errors  of  the  heights  of 
Douglas  fir  used  (right  bar  of  each  pair)  and 
not  used  by  porcupines  in  the  leeward  pure  forest 
on  the  Zoratti  study  area  for  each  of  the  20  years 
prior  to  measurement*  Also  included  for  each  year 
are  the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t-test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Figure  8.  Means  and  standard  errors  of  the  conformation 

indices  of  Douglas  fir  used  (right  bar  of 
each  pair)  and  not  used  by  porcupines  in 
the  leeward  pure  forest  on  the  Zoratti  study 
area  for  each  of  the  20  years  prior  to 
measurement.  Also  included  for  each  year  are 
the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t-test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Figure  9.  Means  and  standard  errors  of  the  diameters  at 

breast  height  of  Douglas  fir  used  (right  bar 
of  each  pair)  and  not  used  by  porcupines  in 
the  leeward  pure  forest  on  the  Skelton  study 
area  for  each  of  the  20  years  prior  to 
measurement.  Also  included  for  each  year  are 
the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t-test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Figure  10.  Means  and  standard  errors  of  the  heights  of 

Douglas  fir  used  (right  bar  of  each  pair)  and 
not  used  by  porcupines  in  the  leeward  pure  forest 
on  the  Skelton  study  area  for  each  of  the  20  years 
prior  to  measurement.  Also  included  for  each  year 
are  the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t-test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Figure  11.  Means  and  standard  errors  of  the  conformation 

indices  of  Douglas  fir  used  (right  bar  of 
each  pair)  and  not  used  by  porcupines  in 
the  leeward  pure  forest  on  the  Skelton  study 
area  for  each  of  the  20  years  prior  to 
measurement.  Also  included  for  each  year  are 
the  number  of  used  and  unused  trees.  An 
asterisk  denotes  a  statistically  significant 
difference  (p  <  0.05)  between  means  as  determined 
by  a  t^test.  All  t-tests  employed  a  pooled 
variance  estimate. 
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Table  23.  Regression  statistics  describing  attempts  to 

explain  year  to  year  variation  in  the  average  size 
of  Douglas  fir  used  by  porcupines  in  the  leeward 
pure  forests  on  the  Zoratti  and  Skelton  study 
areas  (see  Figures  6-1 1j  ,  Independent  variables 
that  could  have  entered  the  equation  included 
use  age,  the  number  of  plots  used  per  year  and 
the  number  of  trees  used  per  year.  Additional 
statistics  have  been  provided  in  Appendix  3. 


Independent 
Variable^1*  Beta 

R2 

P 

Weighted 
Mean±se<  2  J 

Zoratti 

in=17) 

DBH 

Use  Age  {3}  -0.831 

0.  690 

p<0.001 

1 6 • 1 1 ±0 .  645  cm 

Ht 

p>0 . 05 

9 . 12±0 . 298  m 

CI  {1} 

Use  Age  {4}  -0.487 

0.238 

0 . 05>p>0 . 025 

1 . 75±0  .056 

Skelton 

(n=16) 

DBH 

p>0 • 05 

25. 96±0. 799  cm 

Ht 

p>0 . 10 

12 . 00±0 . 40 8  m 

CI  {4} 

Use  Age  {4}  0.526 

0.277 

0  •  05>p>0 . 025 

2 . 24±0 . 056 

1  Numbers  in  braces  refer  to  the  type  of  transformation 
applied:  {1}  squared;  {2}  square  root;  {3}  log;  {4}  inverse. 

2If  age  explained  a  significant  proportion  of  the  variation 
in  the  particular  attribute,  the  corresponding  mean  was 
adjusted  so  that  the  effect  of  age  between  years  is 
negligible. 
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Trees  used  by  porcupines  in  the  Skelton  leeward  forest, 
during  the  20  years  prior  to  measurement,  were  generally 
larger  than  the  trees  used  on  the  Zoratti  area  during  the 
same  period  (DBH  t=9.647  p<0.001  P  df=31:  height  t=5. 750 
p<0.001  P  df =3 1 :  Cl  t=6. 156  p<0.001  P  df=31:  Table  23).  If 
porcupines  exercised  size  preferences  solely  related  to  the 
ease  of  climbing  trees,  local  differences  in  the  size  of 
used  trees  of  the  magnitude  observed  would  not  be  expected 
unless  the  preferred  size  was  unavailable.  Although  trees 
on  the  Skelton  area  were  generally  larger  than  on  the 
Zoratti  area  (Table  8) ,  this  difference  is  not  sufficient  to 
explain  the  disparity  in  the  size  of  used  trees.  The 
relatively  greater  intensity  of  porcupine  feeding  on  the 
Zoratti  area  (Table  14)  may  have  forced  some  animals  to  use 
trees  smaller  than  the  preferred  size.  Storm  and 
Halvorson  (1967)  presented  convincing  evidence  that  the 
removal  of  bark  by  porcupines  greatly  affects  the  vigour  of 
a  tree  for  up  to  ten  years  so  that  in  a  forest  that  has  been 
heavily  used  for  several  years,  the  number  of  rapidly 
growing  trees  within  the  preferred  size  class  would  be 
reduced  and  porcupines  would  be  faced  with  a  choice  between 
slower  growing  trees  of  the  appropriate  size  or  smaller 
vigorous  trees.  Election  between  these  two  alternatives 
would  depend  upon  the  relative  significance  of  tree  size  and 
vigour  to  porcupines  and  upon  a  porcupine's  ability  to 
recognize  indicators  of  vigour  such  as  the  amount  of  annual 
branch  extension  or  the  quantity  of  recent  needles. 
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Temporal  Changes  in  the  Intensity  of  Use 

The  number  of  plots  used  per  year  can  be  regarded  as  an 
index  of  the  number  of  porcupines  present  on  a  study  area 
per  winter.  Considering  this  index  of  porcupine  numbers  for 
a  period  of  several  decades  should  provide  an  indication  of 
the  number  of  porcupines  inhabiting  the  study  areas  through 
time.  This  index  is  probably  more  reliable  than  an  index 
based  on  the  number  of  individual  trees  used  per  year  since 
individual  porcupines  tend  to  concentrate  their  feeding* 
during  a  single  winter*  in  a  very  small  area  (Spencer  1964) 
and  the  number  of  used  trees  counted  for  a  particular  winter 
would  be  greatly  affected  by  the  location  of  the  sampling 
plots  with  respect  to  these  "feeding  pockets". 

Figure  12a  illustrates  that  the  intensity  of  use  of  an 
area  by  porcupines  is  not  consistent  through  time.  If  the 
short-term  fluctuations  are  smoothed  by  plotting  three  year 
means  and  the  curves  are  shifted  so  that  the  mean  ages  of 
the  forests  are  equivalent*  the  two  curves  follow  one 
another  quite  closely  (Figure  12b:  r=0*922  r2=0. 850  p<0.001 
n=14) .  The  final  drop  of  the  Skelton  curve  (Figure  12b) 
would  seem  to  predict  a  similar  decline  in  porcupine  numbers 
on  the  Zoratti  study  area.  During  the  winter  (1974-1975) 
following  the  collection  of  the  plot  data  from  which  the 
Zoratti  curve  was  drawn  only  three  porcupines  were  found. 

No  porcupines  could  be  located  on  this  area  during  the 
subsequent  winter  (1975-1976).  This  decrease  in  porcupine 
numbers  may*  however*  have  been  premature  because  the 
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Figure  12.  The  relationship  between  time  and  the  intensity 

of  feeding  by  porcupines  on  the  Zoratti  and 
Skelton  study  areas.  The  upper  graph  (a)  depicts 
the  number  of  plots  used  on  each  area  per 
year  and  the  lower  graph  (b)  illustrates  the 
numbers  of  plots  used  in  relation  to  the  mean 
age  of  the  forest. 


NO.  OF  PLOTS  NO.  OF  PLOTS 
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landowner  of  this  study  area  had  been  shooting  porcupines 
for  several  years  prior  to  the  winter  of  1974-1975, 

A  similar  trend  is  also  evident  in  the  data  presented 
by  Crouch  (1969)  and  Krefting  et  al.  (1962),  although  the 
timing  of  the  peak  in  use  varies.  Spencer  (1964)  described 
three  peaks  in  porcupine  use,  separated  by  20  to  30  years, 
from  a  large  area  of  pinon  pine  cPinus  edulis)  forest  in 
south-western  Colorado,  but  did  not  reveal  whether  these 
peaks  were  experienced  in  all  parts  of  the  forest  or  in 
stands  of  a  similar  size  or  age. 

Changes  in  the  average  number  of  porcupines  feeding  on 
both  study  areas  could  have  been  completely  independent  of 
the  quality  of  the  habitat.  Nothing  is  known  of  the 
activities  of  these  porcupines  during  any  season  but  the 
winter,  although  numerous  studies  have  reported  that 
porcupines  frequent  different  areas  during  summer  and  winter 
(Gabrielson  1928;  Taylor  1935;  Curtis  and  Wilson  1953; 
Brander  1973) .  The  condition  of  summer  habitats  could 
therefore  affect  the  number  of  porcupines  in  winter  habitats 
(see  Fretwell  1972),  but  the  close  association  between  the 
number  of  plots  used  per  year  and  the  mean  ages  of  the 
forests  (Figure  12b)  suggests  that  the  intensity  of  use  of 
the  winter  habitats  was  dependent  on  their  quality. 

These  observations  support  an  earlier  suggestion  that 
the  intensity  of  use  in  a  forest  should  rise  and  fall  as  the 
size  of  vigorous  trees  passed  through  the  range  of  sizes 
preferred  by  porcupines.  The  observed  decline  in  the 


•  •  ,A 

» 

r 


' 


90 


intensity  of  porcupine  use  of  both  study  areas  could  have 
further  resulted  from  the  concurrent  reduction  of  vigour 
associated  with  ageing  in  the  entire  forest.  An  additional 
implication  of  these  findings  is  that  local  changes  in  the 
basic  population  size  of  porcupines  may  be  the  result  of 
selection  of  a  "better”  habitat  rather  than  changes  in 
natality  or  mortality. 
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Characteristics  of  Porcupine  Feeding 

Repeated  Feeding  in  Individual  Trees 

Use  of  a  particular  tree  during  a  given  winter  may 
influence  its  further  use  by  porcupines  in  subsequent  years. 
Several  studies  have  reported  that  porcupines  return  to  the 
same  tree  year  after  year,  often  enlarging  old  scars 
(Curtis  1941;  Shapiro  1949;  Spencer  1964).  Although 
porcupines  may  select  trees  that  have  been  used  in  the  past, 
repeated  use  of  individual  trees  could  have  been  a  random 
event.  The  latter  suggestion  can  be  tested  by  examining  the 
goodness  of  fit  between  a  Poisson  distribution  and  the  group 
frequencies  of  trees  classed  according  to  the  number  of 
years  that  they  were  fed  upon  by  porcupines  (see  Sokal  and 
Rohlf  1969). 

With  the  exceptions  of  Douglas  fir  in  the  leeward  pure 
forests,  the  trees  on  both  study  areas  were  used  from  year 
to  year  in  a  random  fashion  (Tables  24  and  25) .  This 
implies  that  even  though  there  was  a  preferred  type  of  tree, 
as  shown  earlier  (Tables  18,  19  and  20),  there  were  enough 
of  that  type  that  they  could  be  used  in  a  manner  independent 
of  previous  use.  Curtis  and  Wilson  (1953)  also  examined 
this  relationship  statistically  and  came  to  a  similar 
conclusion. 

Douglas  fir  in  the  leeward  pure  forests  were  used  in 
contrasting  patterns  on  the  two  study  areas.  Individual 
trees  in  low  density  stands  on  both  areas  were  used  during 
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Table  24.  The  frequency  of  repeated  use  of  individual  trees 


on  the  Zoratti  study  area.  The  log  likelihood 
ratio  tests  (G)  are  based  on  a  comparison  between 
the  observed  frequencies  if)  and  the  expected 
frequencies  (f)  of  a  similar  population  (i.e. 
identical  mean)  with  a  Poisson  distribution. 

The  degrees  of  freedom  associated  with  this 
goodness  of  fit  test  are  the  number  of 
classes  -  2  iSokal  and  Rohlf  1969),  so  that  this 
test  is  inappropriate  for  less  than  three  classes. 


93 


Q> 

1 

1 

in 

M 

» 

CM 

1 

14 

o 

3 

» 

• 

• 

• 

• 

• 

1 

OV 

CM 

o 

04 

1 

CM 

to 

■d- 

CN 

CM 

1 

<T> 

O 

• 

M<4H 

» 

CM 

lO 

lO 

r— 

» 

vO 

VO 

o 

*0  'H 

1 

CM 

r— 

1 

» 

• 

V 

U  &4 

• 

1 

00 

o 

vo 

CM 

04 

<tJ 

1 

1 

r— 

* 

» 

CM 

5J- 

Ov 

f— 

CM 

t 

^9- 

<1)  44 

1 

O 

CO 

r- 

1 

<D 

♦ 

CM 

T— 

1 

1 

f 

» 

1 

« 

» 

r- 

in 

I 

1 

TJ 

» 

• 

• 

• 

t 

o> 

in 

r— 

U  rH 

1 

r- 

CD 

vO 

1 

CO 

co 

• 

(t)  <d<4-4 

1 

vo 

VO 

1 

T— 

m 

O 

*  4-> 

1 

CO 

t 

• 

• 

A 

*o  o 

1 

» 

r-~ 

o 

T— 

T— 

04 

a  eh 

» 

1 

CO 

*H 

1 

CO 

in 

1 

ts  4-1 

1 

to 

t"* 

1 

1 

1 

co 

1 

1 

ft 

f 

1 

1 

* 

vO 

vO 

f 

» 

1 

*d 

I 

• 

• 

• 

1 

CM 

r» 

T~ 

I 

to 

co 

to 

1 

OV 

CM 

• 

(d  0><44 

1 

o 

ID 

» 

r— 

vO 

o 

*  C 

I 

CO 

» 

• 

• 

A 

h 

1 

t 

O 

O 

r~ 

r— 

O. 

C  04 

1 

1 

r* 

■iH 

1 

CO 

CO 

I 

CO 

t*  4-1 

1 

o 

vO 

1 

1 

ft 

CO 

• 

ft 

f 

f 

ft 

f 

? 

ft 

I 

» 

o 

o 

<Tt 

f 

1 

no 

* 

• 

• 

• 

1 

CT> 

(4 

♦ 

VO 

o 

o 

» 

<t)  <4-1 

I 

to 

t— 

1 

r— 

»  M 

» 

1 

• 

no  <h 

1 

1 

r* 

O 

*4 

0  t*4 

1 

1 

vO 

•H 

1 

vO 

O 

f— 

• 

tx  4-1 

t 

1 

» 

ft 

to 

T" 

1 

t 

1 

ft 

W 

f 

l 

l 

1 

f 

I 

1 

1 

H  nO 

l 

• 

#  m  d) 

l 

o 

r— 

CM 

ro 

» 

c 

IX 

U 

44 

04 

<D  W 

t 

1 

*0 

>H  CD 

l 

I 

(continued) 


■ 


. 


I 


i 

O' 

CO 

04 

04 

CO 

f 

N 

in 

i  d 

O 

3 

t 

• 

• 

• 

• 

■ 

4 

O 

r— 

o 

i  > 

u 

(0 

i 

r— 

CD 

CO 

T” 

r~ 

1 

vO 

CO 

• 

4 

<y 

»H 

»“K3 

i 

O'. 

« — 

O' 

T~ 

» 

3 

O 

o 

4  3 

N 

3  » 

<d 

» 

VO 

no 

» 

• 

• 

A 

4  cd 

d 

3 

i 

* 

vO 

O 

LO 

OJ 

Oi 

•  3 

Q> 

3  O'* 

o 

i 

1 

O' 

1  u 

> 

(0  vO 

Eh 

t 

3- 

04 

C" 

Ol 

4 

o 

4  CD 

d 

av 

3 

i 

CT> 

OJ 

vO 

«“ 

» 

r— 

4  a. 

<D 

CO  vr- 

i 

VO 

CO 

t 

4  X 

•H 

» 

1 

1  0) 

3  3 

• 

1 

4 

'd 

3  *H 

i 

1 

4  a 

•H 

3 

3 

i 

1 

1  d 

3 

Q)  O 

0) 

» 

T~ 

00 

3 

» 

N 

4 

O 

3  « 

X 

i 

• 

• 

• 

• 

♦ 

Oi 

T“ 

4  3 

» 

d 

•H 

» 

r- 

CO 

vO 

o 

1 

CO 

CO 

• 

4  d 

3  <3 

se 

» 

T— 

CO 

1 

CO 

m 

o 

4  -3 

w 

3  3 

d 

i 

T— 

• 

• 

• 

A 

4 

w 

O  d 

3 

3 

i 

» 

CO 

o 

OI 

*— 

Cu 

4  CO 

d 

*3 

u 

o 

i 

1 

vO 

4  CO 

d  r- 1 

d 

E-t 

i 

r- 

04 

00 

04 

1 

T— 

4  d 

U 

u  d 

» 

3 

♦ 

04 

CO 

* 

4  3 

M 

cu 

i 

r- 

» 

4  U 

■P 

O  O 

<d 

r 

1 

4 

3 

3  CO 

3 

i 

1 

4  0) 

<D 

— * 

t 

4 

1  3 

U 

3 

i 

1 

4  O 

d 

cu  a) 

3 

i 

1 

• 

TO 

3  > 

<D 

i 

m 

CO 

cr> 

CO 

4 

4  3 

»d 

X 

» 

• 

• 

• 

1 

3 

4  CO 

d 

3  3 

>H 

• 

r- 

o 

04 

o 

» 

00 

4  d 

a 

S3 

«4-l 

i 

f» 

o» 

1 

OJ 

4  a) 

<D 

O  3 

CD 

i 

1 

• 

4  H 

3 

O  0) 

3 

a 

» 

1 

if) 

© 

«4 

4 

•P 

0) 

(-i 

•H 

« 

» 

CO 

4  3 

3  U 

d 

Cu 

i 

VO 

CO 

3 

04 

1 

i  d 

3 

0)  X 

s 

3 

i 

r- 

I 

4 

P> 

a)  a> 

CD 

« 

4 

4  0) 

3 

CD 

i 

1 

i  in 

3 

PI 

i 

1 

4  3 

<u  3 

• 

1 

4  d 

W 

>  3 

t 

4 

4  U 

W 

d  O 

3 

♦ 

4 

4  V 

d 

3  » 

<D 

» 

00 

o 

VO 

4 

4  3 

X 

i 

• 

• 

• 

4 

CO 

4 

U 

CO  >i 

•H 

i 

3 

00 

co 

4 

o 

4  3 

0)  u 

E 

<3 

i 

3 

<r" 

4 

3 

4  0) 

W 

CO  3 

M 

i 

4 

• 

4  E 

•H 

co  <y 

3 

•H 

t 

4 

O' 

o 

4  U 

3 

d  s 

u 

Cm 

t 

4 

vO 

4  O 

3 

r— i  tr* 

d 

i 

3 

<7i 

3 

4 

4  3 

o  a> 

» 

3 

i 

3 

r— 

4 

4  3 

dv 

M 

CD 

i 

4 

4  CD 

3 

3  3 

CD 

i 

4 

4  Cu 

•H 

3 

3 

» 

4 

4 

CU 

a)  3 

i 

4 

4  3 

e 

• 

3  0) 

3 

i 

4 

4  O 

3 

s 

d*3 

0> 

CD 

i 

co 

3* 

av 

vO 

4 

(VI 

1  3 

•H 

o 

a>  u 

3 

3 

i 

• 

• 

• 

• 

4 

OJ 

3 

* — 

4 

'd 

M  Q 

a 

3 

« 

vO 

o* 

a> 

co 

4 

o 

OJ 

• 

1  W 

d 

0)<U  Oi 

•H 

a, 

cd<3 

i 

T" 

T— 

4 

T“* 

r— 

o 

4  d 

3 

<u 

X 

3 

o 

i 

4 

• 

• 

A 

•  a 

d 

4-4 

3  9 

C 

3 

3 

i 

4 

r- 

3 

r— 

a. 

r 

3 

CO 

o 

(-1 

(-4 

» 

4 

3 

4  3 

in 

d  CD 

u 

d 

cu 

» 

CO 

04 

o 

3 

4 

4  CO 

3 

a)  3 

» 

to  3 

i 

T— 

04 

r— 

4 

4  CP 

3 

O  r|  3 

« 

0) 

i 

4 

4  -M 

d 

3- 

CD 

i 

4 

4  4 

3 

W 

O  3 

04 

•-4 

» 

4 

4  C5 

3 

<y 

»  d 

i 

4 

4 

Q> 

3  3 

<1> 

co 

» 

4 

4  OJ 

CO 

3 

3 

<— 4 

M  3 

i 

4 

4  ’H 

CO 

C7V  Q) 

-Q 

=** 

d  a) 

i 

o 

04 

CO 

3 

4 

3 

\x 

O 

3 

a. 

•  3 

CD 

0)  3  O 

d 

cd  to 

i 

4 

3 

4  H 

<— i  -d 

Eh  W 

E-* 

>H  C3 

i 

4 

4 

N 

• 

- 

' 

95 


> 


’ 

* 

. 

■ 

■ 


' 


Table  25.  The  frequency  of  repeated  use  of  individual  trees 
on  the  Skelton  study  area.  The  log  likelihood 
ratio  tests  (G)  are  based  on  a  comparison  between 
the  observed  frequencies  if)  and  the  expected 
frequencies  (f)  of  a  similar  population  (i.e. 
identical  mean)  with  a  Poisson  distribution. 

The  degrees  of  freedom  associated  with  this 
goodness  of  fit  test  are  the  number  of 
classes  -  2  iSokal  and  Rohlf  1969),  so  that  this 
test  is  inappropriate  for  less  than  three  classes. 
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‘This  G-test  was  not  performed  because  at  least  one  class  had  an  expected 
value  less  than  5  and  lumping  this  class  with  the  adjacent  class  would 
leave  zero  degrees  of  freedom. 

2The  four  least  frequent  classes  have  been  combined  to  calculate  this 
statistic,  so  that  the  expected  frequency  would  exceed  five. 
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more  years  than  trees  in  high  density  stands  (Table  26: 
Zoratti  G=  1 0 . 6  08  p<0.001;  Skelton  G=16.360  p<0.001  df=1), 
suggesting  that  the  preferred  type  of  tree  was  uncommon 
relative  to  the  number  of  porcupines.  Because  of  this 
relative  scarcity,  desirable  Douglas  fir  on  the  Skelton  area 
had  a  high  probability  of  being  used  during  several  years 
(Table  25)  and  trees  in  low  density  stands  were  used,  during 
at  least  one  year,  more  frequently  than  trees  in  high 
density  stands  (G=104. 676  p<0.001  df =1 ) .  On  the  Zoratti 
area  however,  the  number  of  resident  porcupines  was 
apparently  so  great,  during  many  winters,  that  trees  in  high 
and  low  density  stands  were  used,  at  least  once,  in 
equivalent  proportions  (G=2.8981  p>0.05  df=1)  .  If  used 
trees  in  high  density  stands  on  the  Zoratti  area  were  of 
"poor  quality",  the  probability  of  their  being  used  during  a 
subsequent  winter  would  have  been  lower  than  for  "good 
quality"  trees,  resulting  in  more  trees  than  expected  being 
used  only  once  (Table  24) . 

Porcupine  Feeding  and  Stand  Density 

The  demonstrated  correlation  between  stand  density  and 
the  proportion  of  trees  used  by  porcupines  (Table  17)  and 
the  relationship  between  age  and  the  intensity  of  use 
(Figure  12b)  suggest  two  possible  patterns  of  forest 
utilization.  When  the  forest  was  relatively  young  and 
frequented  by  a  limited  number  of  porcupines,  low  density 
(preferred)  stands  should  have  received  the  most  attention. 
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Table  26,  A  comparison  of  the  frequency  of  porcupine 
feeding  in  high  and  low  density  plots.  The 
numbers  in  parentheses  are  row  percentages. 


Zoratti 

Number  of  Years  Used 


Density 

(#trees/100  sq.  m)  0  1  >1 


127  95  39 

<11.00  (48.  7)  (36.  4)  (14.9) 

75  99  13 

>11.00  (40.  1)  (52.9)  (7.0) 


Skelton 


Density 

(#trees/100  sq.  m) 

Number 

of  Years 

Used 

0 

1 

>1 

206 

119 

55 

<11.00 

(54. 2) 

(31.3) 

<14.5) 

>11.00 

7  03 

128 

18 

(82.8) 

(15.1) 

(2.1) 

• 

• 

99 


As  the  forest  grew  and  porcupine  numbers  increased,  the 
porcupines  could  have  crowded  into  the  preferred  type  of 
stand,  or  some  animals  could  have  moved  into  the  less 
preferred  stands  (see  Fretwell  1972) . 

Whereas  porcupines  on  the  Zoratti  area  used  a  similar 
proportion  of  trees  in  high  and  low  density  plots  (see 
Table  14:  G=1.586  p>0.1  df=1),  porcupines  on  the  Skelton 
area  showed  a  preference  for  trees  in  low  density  leeward 
plots  (see  Table  14:  G=93.130  p<0.001  df=1) .  During  the 
32  year  period  from  1942  to  1974,  an  average  (±  se)  of 
29.9  ±  3.78  percent  of  the  leeward  pure  plots  on  the  Zoratti 
area  were  used  per  year,  while  on  the  Skelton  area  during  a 
similar  period  (1940-1972),  17.0  ±  1.78  percent  of  the 
leeward  plots  were  used  (t=:2.968  0.  01>p>0.001  S  df=44)  .  The 
high  intensity  of  use  experienced  on  the  Zoratti  area  may 
not  have  allowed  the  expression  of  a  preference  for  stands 
with  less  than  11  trees/100  sg .  m. 

Figure  13  illustrates  an  increase  in  the  use  of  high 
density  plots  with  increasing  porcupine  numbers  on  the 
Skelton  area  and  provides  additional  evidence  that  high 
density  stands  are  poorer  quality  porcupine  habitat.  (Only 
the  period  from  1949  to  1970  (see  Figure  12a)  was  considered 
because  the  chance  use  of  a  high  density  plot  when  the 
population  was  small  would  have  a  much  greater  effect  on  the 
proportion  of  used  plots  that  were  high  density  than  during 
population  peaks.)  In  addition,  this  observation  suggests 
that  individual  porcupines  were  aware  of  the  distribution  of 
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Figure  13.  The  relationship  between  the  number  of  porcupines 

feeding  in  the  Skelton  leeward  forest,  as  measured 
by  the  number  of  plots  used  per  year,  and  the  use 
of  high  density  stands. 


HIGH  DENSITY  PLOTS/TOTAL  PLOTS  USED  PER  YEAR 
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other  porcupines  around  them,  even  though  they  are  generally 
solitary  animals  (Taylor  1935;  Curtis  and  Kozicky  1944). 

Location  and  Amount  of  Feeding  Within  Individual  Trees 

Although  porcupines  feed  on  both  the  inner  bark  and 
foliage  of  conifers,  evidence  of  feeding  on  foliage  is 
ephemeral  and  cannot  be  adequately  assessed  after  several 
years.  Porcupines  tend  to  feed  in  the  upper  portions  of 
conifers  on  both  the  bole  and  lateral  branches  (Tables  27 
and  28a:  Gabrielson  1928;  Taylor  1935;  Shapiro  1949) . 

Feeding  scars  on  the  upper  surfaces  of  branches  were  usually 
encountered  over  a  greater  height  range  in  a  particular  tree 
than  were  scars  on  the  bole. 

Usually  less  than  650  sg.  cm  of  bark  were  removed  from 
an  individual  tree,  although  many  trees  suffered  greater 
damage  (Table  28b) .  Feeding  most  commonly  occurred  on 
portions  of  the  tree  unprotected  by  corky  bark,  such  as  the 
distal  ends  of  branches  and  the  upper  bole.  Porcupines 
often  succeeded  in  removing  bark  from  the  entire 
circumference  of  the  bole  of  a  tree  resulting  in  the  death 
of  the  upper  crown  (Table  28c) .  This  type  of  damage  has  the 
greatest  consequence  on  tree  vigour  (Storm  and 
Halvorson  1967)  and  often  can  alter  the  form  of  an 
individual  tree,  particularly  if  the  tree  is  regirdled  in 
later  years.  Trees  in  low  density  stands  in  the  leeward 
pure  forest  often  had  more  bark  removed  by  porcupines  during 
single  winters  than  trees  in  high  density  stands  (Zoratti 
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Table  27.  The  average  (±  standard  deviation)  position  of 
porcupine  feeding  within  trees  on  the  Zoratti 
and  Skelton  areas.  Except  for  the  column  on 
the  extreme  right,  position  is  measured  on  a 
scale  on  which  ground  level  has  a  value  of  zero 
and  the  top  of  the  tree  has  a  value  of  one. 
Calculation  of  the  position  of  the  mode  within 
the  feeding  range  involved  a  similar  scale  on 
which  the  bottom  and  top  of  the  feeding  range 
within  a  tree  had  values  of  zero  and  one, 
respectively . 


Mode 


Top 

of  Feeding 
Range 

Bottom 
of  Feeding 
Range 

Mode 

of  Feeding 
Within  Tree 

of  Feeding 
Within 

Feeding  Range 

ZWF1 

0.51±0.204 

111)2 

0 . 29±0 .212 
(11) 

0.38±0.207 

(9) 

0.6  2±0 . 243 
(9) 

ZPF 

0. 80±0 .200 
(2  33) 

0 . 37±0 • 203 
(233) 

0.65±0.190 

(121) 

0.59±0. 181 
(165) 

ZMF 

0. 58±0 .275 
1244) 

0. 24±0. 170 
124) 

0.40±0.214 

(16) 

0 . 60±0 . 164 

(13) 

ZPS 

0*78 ±0.131 
143) 

0.  29± 0. 130 

143) 

0. 53±0 . Ill 
(16) 

0. 59±0. 148 
(18) 

ZWP 

0.63±0.275 

149) 

0. 33± 0.242 
(49) 

0. 5 1±0. 23  7 
140) 

0 . 57±0 . 242 

140) 

ZMP 

0. 66±0. 249 

(16) 

0. 42±0. 199 

(16) 

0 . 57±0 . 1 09 
(7) 

0. 53±0. 182 

(12) 

SWF 

0.78+0.176 

(19) 

0. 37±0. 199 

(19) 

0.67±0.171 

(10) 

0.7  2±0 .171 
(19) 

SPF 

0. 82±0 .197 
(179) 

0. 49±0. 222 
1179) 

0. 66±0. 180 
(78) 

0.59±0. 244 
(182) 

SWP 

(2) 

(2) 

(2) 

0 . 67±0 . 3 1 3 

(6) 

Explanation  of  this  comparison  code  is  provided  in  Table  9 
2Sample  size. 
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Table  28.  Some  characteristics  of  porcupine  feeding  in 

individual  trees  on  the  Zoratti  and  Skelton  areas. 
Each  entry  indicates  the  number  of  trees  in  the 
respective  class.  Part  a)  is  based  upon  trees 
that  were  used  during  a  single  year  only. 
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G=22. 8  72  p<0.001  d£=2:  Skelton  G=18.  942  p<0.001  df=2) . 

The  porcupines  on  the  Zoratti  and  Skelton  areas  appear 
to  have  fed  on  the  portions  of  trees  in  which  food  was  most 
easily  attainable  and  probably  most  nutritious.  The  phloem 
and  foliage  of  conifers  must  supply  much  of  the  nutrients 
required  for  spring  growth  and  therefore  contain  the  highest 
concentrations  of  fats  and  carbohydrates  of  the  aerial 
tissues  (see  Kramer  and  Kozlowski  1960  for  a  review).  Tree 
growth  is  annually  initiated  at  the  bud  bases  and  nutrients 
are  drawn  from  the  adjacent  tissues  (Kozlowski  and 
Keller  1966)  so  that  the  bark  and  branches  within  the  crown 
are  important  storage  structures.  The  tendency  of 
porcupines  to  feed  on  foliage  from  young  branches 
(Taylor  1935;  Curtis  1941;  Shapiro  1949)  and  inner  bark 
within  the  tree  crown  probably  resulted  in  the  consumption 
of  the  most  nutritious  food  available. 
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CONCLODING  DISCUSSION 

Preferences  and  Selection 

Whether  porcupines  actually  selected  their  food  and 
habitat  in  the  Zoratti  and  Skelton  forests  or  merely  relied 
upon  chance,  is  central  to  understanding  the  patterns  of 
forest  use  which  have  been  observed.  The  disproportionate 
utilization  of  a  specific  segment  of  a  resource  relative  to 
its  occurrence  in  the  environment  has  been  regarded  as  a 
preference  throughout  the  preceding  discussion.  Selection, 
the  active  expression  of  preferences  involving 
discrimination  between  alternatives,  is  commonly  inferred 
from  a  demonstration  of  preferences,  often  erroneously 
iWiens  1976) .  The  apparent  relationship  between  the 
location  of  a  forest  and  the  intensity  of  its  utilization  by 
porcupines  is  a  germane  example  of  the  occurrence  of 
preferences  without  the  involvement  of  selection. 

Changes  in  the  preferences  of  porcupines  in  relation  to 
the  number  of  individuals  using  an  area  intimate  selection. 
The  increased  feeding  in  high  density  stands  on  the  Skelton 
area  when  porcupines  were  relatively  numerous  iFigure  13) 
implies  that  these  stands  were  accessible  and  usable,  but 
probably  not  as  desirable  as  low  density  stands.  Deviation 
from  randomness  in  the  repeated  use  of  individual  trees  in 
the  leeward  pure  forests  on  both  study  areas,  even  though 
feeding  in  trees  in  other  communities  was  a  random  event 
(Tables  24  and  25) ,  is  also  suggestive  of  selective 
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behaviour.  These  observed  preferences  were  seemingly 
mediated  by  the  porcupines  themselves,  indicating  an  ability 
on  their  part  to  detect  differences  in  quality. 

Consumers  cannot  afford  to  be  selective  when  resources 
are  limited  (Emlen  1966,  1968)  and  the  presence  of 
selection- based  preferences  therefore  implies  a  relative 
abundance  of  food  and  habitat.  Relaxation  of  selectivity 
was  probably  necessary  on  the  more  extensively  used  Zoratti 
area  resulting  in  the  apparent  preference  for  smaller  trees 
than  were  used  on  the  Skelton  area.  The  characteristics  of 
the  trees  used  by  porcupines  on  the  Skelton  area  may 
consequently  be  a  better  reflection  of  the  physical 
qualities  considered  by  porcupines  in  selecting  trees. 

Food  Habits  and  Habitat  Selection 

Winter  food  and  habitat  selection  by  porcupines  are 
apparently  strongly  interdependent.  If  a  patch  of  bark  or  a 
branch  of  foliage  are  considered  the  food  items  of  a 
porcupine  feeding  within  a  tree,  the  tree  itself  assumes  the 
characteristics  of  a  repeatable  but  uneven  microhabitat 
throughout  which  the  food  items  are  scattered.  Each 
individual  tree  supports  many  potential  food  items,  although 
some  are  probably  nutritionally  or  energetically 
satisfactory  but  unattainable  because  of  their  position.  A 
porcupine  that  is  foraging  efficiently  should  therefore 
select  trees  on  the  probability  of  obtaining  a  relatively 
large  quantity  of  suitable  food  items  (Royama  1970; 
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MacArthur  1972)  . 

This  type  of  selection  would  require  the  presence  of 
characters,  such  as  bole  diameter  or  crown  width,  which 
could  be  used  to  reliably  distinguish  profitable  trees, 
preferably  without  unnecessarily  expending  energy  in 
climbing.  Chemical  characteristics  may  be  of  little 
assistance  in  the  selection  of  trees  from  ground  level, 
because  the  living  tissues  at  the  base  of  a  tree  are  often 
protected  by  several  centimeters  of  corky  bark  which  could 
hinder  assessment  of  the  present  nutrient  quality  of  the 
tree.  With  the  exception  of  indications  of  porcupines 
feeding  at  ground  level  on  limber  pine  in  windward 
communities  (see  also  Gill  and  Cordes  1972),  basal  scars 
were  infrequently  encountered  during  this  study. 

In  addition  to  food,  habitat  must  provide  protection 
from  predators  and  the  elements.  Rock  dens  supply  these 
requirements  when  they  are  available  to  porcupines 
(Curtis  1941;  Shapiro  1949),  but  in  their  absence,  trees, 
thickets  of  brush,  and  windfalls  serve  as  shelter 
(Taylor  1935;  Curtis  and  Wilson  1953).  Because  neither 
study  area  offered  suitable  rock  dens,  individual  trees  must 
have  furnished  both  food  and  shelter  for  the  resident 
porcupines . 

The  selection  of  individual  trees  by  porcupines  could 
have  occurred  by  two  alternate  procedures.  A  porcupine 
could  seek  out  the  appropriate  community  and,  having  found 
it,  search  for  a  suitable  stand  and  finally  decide  upon  an 
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attractive  tree;  or  the  entire  selection  process  could 
involve  only  the  search  for  an  acceptable  tree.  Although 
either  method  of  selection  could  have  resulted  in  the 
observed  preferences  because  large,  vigorous,  open-grown 
trees  were  most  abundant  in  low  density  stands  in  the 
leeward  pure  forests,  the  occurrence  of  inter-community 
preferences  and  the  absence  of  interspecific  preferences 
suggest  porcupines  use  the  more  efficient  procedure  of 
examining  successively  smaller  units. 

The  strategy  of  food  selection  adopted  by  a  particular 
individual  or  species  is  determined  by  the  interaction  of 
various  constraints  imposed  by  time,  energy  and  nutrient 
requirements  (Ellis  et  al.  1976).  Porcupines  are  probably 
not  confronted  with  a  need  to  minimize  the  time  spent 
feeding  during  the  winter  because  they  enjoy  protection  from 
most  predators  while  situated  in  trees  (Taylor  1935)  and 
they  are  not  involved  in  other  time  consuming  activities. 

In  addition,  satisfaction  of  energy  requirements  is 
presumably  not  difficult  for  porcupines  because  these 
animals  are  often  inactive  for  long  periods  during  the 
winter  (Taylor  1935;  Brander  1973)  and  they  can  withstand 
moderately  cold  temperatures  without  dramatically  elevating 
their  metabolic  rate  (Irving  et  al.  1955:  Winter  lower 
critical  temperature  of  Alaskan  porcupines  was  -12°  C) .  The 
fulfillment  of  nutrient  requirements,  particularly  in 
forests  that  offer  a  limited  number  of  food  species,  may 
therefore  be  the  most  influential  factor  in  determining  the 
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selection  of  winter  food  by  porcupines. 

Food  value  is  the  ultimate  determinant  of  selection 
behaviour,  even  if  the  characters  considered  during  the 
selection  process  are  only  correlated  with  a  food's 
nutritive  qualities  (Fretwell  1972) .  The  quality  of  a  food 
item  is  a  complex  property  which  involves  not  only  its 
nutrient  content,  but  also  the  digestion  and  assimilation 
efficiency  characteristics  of  the  consuming  individual 
(Longhurst  et  al.  1968)  which  are  subject  to  the  influence 
of  secondary  plant  compounds  (Freeland  and  Janzen  1974).  In 
addition,  an  animal*  s  response  to  a  particular  stimulus, 
such  as  a  food  with  a  specific  nutrient  content,  can  be 
greatly  modified  by  its  physiological  state  (Cabanac  1971). 
An  awareness  by  porcupines  of  chemical  qualities  is  not  an 
assumption  required  to  explain  the  observed  porcupine 
preferences  of  communities,  stands  or  individual  trees, 
because  vigorous  trees  can  be  identified  by  obvious  physical 
features.  Chemical  attributes  may,  however,  be  more 
important  within  individual  trees  because  porcupines  appear 
to  select  food  items  through  olfaction  (Murie  1926; 

Taylor  1935)  . 

Changes  in  food  quality  associated  with  ageing  of  the 
forests  (see  Cowan  et  al.  1950)  have  been  suggested  as  a 
cause  of  the  decline  in  use  of  the  Zoratti  and  Skelton  areas 
by  porcupines.  Temporal  variability  in  food  quality  places 
a  premium  on  dispersal  behaviour  because  dispersal  allows 
the  occupation  of  habitats  that  have  recently  become 
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favourable  (Wiens  1976),  Porcupines  are  particularly  mobile 
during  the  spring  and  fall  when  the  availability  of 
herbaceous  food  changes  (Gabrielson  and  Horn  1930; 

Taylor  1935;  Curtis  and  Wilson  1953)  so  that  contact  with 
new  habitats  could  be  frequent.  Whether  the  declines  in 
utilization  of  the  Zoratti  and  Skelton  areas  were 
accompanied  by  a  concurrent  increase  in  the  use  of  adjacent 
forests  is  difficult  to  assess.  Fifteen  to  20  forests 
within  30  kilometers  of  these  study  areas  were  superficially 
examined  for  indications  of  recent  porcupine  activity, 
without  success.  These  forests  were  all  approximately  the 
same  age,  probably  because  of  widespread  fires 
(Dawson  1886) ,  so  that  the  general  quality  of  coniferous 
forests  that  porcupines  use  as  winter  habitat  in  this  region 
could  be  steadily  declining. 

Porcupines  and  Foraging  Theory 

There  exists  in  the  ecological  literature  a  large  body 
of  information  devoted  to  the  optimal  use  of  food  and 
habitat  (see  Ellis  et  al.  1976  and  Wiens  1976  for  reviews) . 
Although  the  theories  proposed  by  individual  studies  often 
include  components  which  are  specific  to  a  limited  type  of 
consumer  (eg.  Schoener  1969;  Krebs  et  al.  1974; 

Pulliam  1974),  basic  precepts  are  apparent  and  they  form  the 
following  allegorical  summary  of  the  observed  dynamics  of 
forest  utilization  by  porcupines. 
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Animals  foraging  optimally  should  specialize  when 
confronted  by  a  relatively  predictable  abundance  of  food 
(Emlen  1966,  1968).  On  a  short-term  basis,  this 
specialization  is  best  accomplished  through  increased 
selectivity.  Abundant  food  often  results  in  a  numerical 
response  on  the  part  of  the  consumer  (Holling  1961),  whether 
through  immigration  or  increased  success  in  reproduction, 
and  an  increase  in  consumer  density.  Continued  exploitation 
of  the  food  supply  and  an  increase  in  the  frequency  of 
behavioural  interactions  with  other  consumers  concurrent 
with  this  increase  in  density  (Holling  1961)  consequently 
reduce  the  relative  suitability  of  the  habitat.  As  this 
decline  in  suitability  continues,  the  consumers  must  become 
less  selective  and  begin  using  previously  less  attractive 
foods  and  habitat  (Fretwell  1972).  The  increased  use  of 
different  habitat  types  is  accompanied  by  a  change  in  the 
distribution  of  consumers  from  aggregation  through 
randomness  to  uniformity  (Grant  and  Morris  1971).  A  decline 
in  the  overall  quality  of  the  habitat,  because  of  over- 
exploitation  or  intrinsic  factors,  should  be  accompanied  by 
a  decline  in  consumer  density. 
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Appendix  1.  Error  in  the  Calculation  of  Past  Diameters, 

The  calculation  of  past  diameters  is  based  upon  the 
equation, 

IRx=  (a+  (b«DBH)  )  •£  (TC-  2  In)  /TC  J  (1) 

in  which:  x  represents  the  number  of  years  prior  to  the 
measurement  of  the  diameter  at  breast  height  (DBH) ;  IR 
represents  the  radius  inside  the  bark  at  time  x;  TC 
represents  the  total  length  of  the  increment  core  at  breast 
height;  I  represents  the  size  of  the  annual  increment 
accrued  during  year  n;  and  a  and  b  are  the  constants  from  a 
linear  regression  equation  relating  the  measured  DBH  and 
inside  radius.  All  quantities  on  the  right  of  this  equation 
were  measured  directly,  with  the  exception  of  TC  which  is  a 
derived  variable  (see  Methods) •  Since  increment  cores  were 
extracted  at  stump  height,  the  core  length  can  be  subdivided 
into  the  equivalent  core  length  at  breast  height  (TC)  and 
the  core  length  differential  (CD)  which  represents  diameter 
growth  before  the  tree  reached  breast  height.  The  age  at 
breast  height  (ABH)  has  been  calculated  using  regression 
techniques  (see  Methods)  and  the  measured  size  of  the  growth 
increment  corresponding  to  this  age  represents  CD.  Error  in 
calculating  past  diameters  is  therefore  primarily  associated 
with  the  estimate  of  ABH  and  perpetuated  through  the 
measurement  of  CD  and  the  calculation  of  TC. 

The  following  examples  are  an  examination  of  the 
magnitude  of  the  error  involved  in  determining  the  past 
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diameters  of  cored  trees  and  the  relationship  of  this  error 
to  tree  age  and  stand  density.  Three  Douglas  fir  trees  from 
the  Zoratti  leeward  pure  association  have  been  used  in  this 
demonstration.  Tree  2304  is  from  a  low  density  plot  while 
3727  and  3749  are  from  the  same  high  density  plot.  In 
addition  to  the  calculated  value  of  TC,  two  values  are 
provided  which  represent  ±  0.5  cm  errors  in  the  estimation 
of  CD. 

Since  TC  contributes  to  both  the  numerator  and  the 

denominator  of  the  quantity  within  the  square  brackets  in 

/ 

equation  (1),  any  error  involved  in  calculating  CD  tends  to 
cancel  itself  out  and  the  influence  of  a  relatively  large 
error  in  the  value  of  TC  on  the  overall  calculation  of  past 
diameters  is  therefore  slight.  The  effect  of  an 
overestimate  is  generally  greater  than  that  of  an 
underestimate.  In  addition,  the  size  of  the  error  is 
inversely  related  to  both  the  age  of  the  tree  and  stand 
density,  but  directly  related  to  the  number  of  years  prior 
to  measurement  of  the  tree  for  which  a  diameter  is  desired. 
Although  no  conclusive  evidence  is  available  it  is  proposed 
that  the  mean  error  associated  with  the  past  diameters  at 
breast  height  does  not  exceed  five  percent  and  that  this 
error  does  not  significantly  affect  the  conclusions 
regarding  the  size  and  form  of  trees  used  by  porcupines. 


. 


Tree  2304 


DBH=28 . 9  cm 
Bark=1.02  cm 
Age=69  yr 
Core=13.22  cm 

1 TC=13. 22-3. 350 

=9.87  cm  10. 

I R  =  (-0. 160+ (0. 456*28. 9) 
=12.0  cm  12. 

DBH  =26.6  cm  27. 

%  Error  0. 

I R  =  (-0. 160+ (0.456*28.9) 
=10.08  cm  10. 

DBH  =22.4  cm  22. 

%  Error  1. 

IR  =(-0. 160+ (0. 456*28. 9) 
=7.40  cm  7. 

DBH  =16.6  cm  17. 


Age  at  Breast  Height=20  yr 
Core  Diff erential=3. 350  cm 

Density=7.77  trees/100  sg. 

cm  9.37  cm 

ll9.87-0.776) /9. 87) 
cm  11.94  cm 

cm  26.5  cm 

0.38 

( (9.87-2 .228) /9.87) 
cm  9.92  cm 

cm  22.1  cm 

1.36 

(9.87-4. 263)/9. 87) 
cm  7.10  cm 

cm  15.9  cm 


37 

)• 

04 

7 

36 

)• 

22 

7 

32 

)  • 

67 

2 


%  Error 


3.  49 


4.40 
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Tree  3727 


DBH=13.3  cm 
Bark=0.76  cm 
Age=73  yr 
Core=5.32  cm 

TC=5.32-2.  735 
=2.59  cm 

IR  =1-0. 160+ (0.456*13. 3) 
=5.71  cm 

DBH  =12.9  cm  1 


%  Error 

IR  =  (-0.  160+ (0.456*13. 3) 
=5.52  cm 


DBH 

=  12.5 

cm 

1 

%  Error 

IR 

=  (-0. 

160+  (0. 

456*13.3) 

=  4.95 

cm 

DBH 

=  11.2 

cm 

1 

Age  at  Breast  Height=20  yr 
Core  Dif ferential=2. 735  cm 

Density=18. 66  trees/100  sq. 

.09  cm  2.09  cm 

•  ((2.59*0.083)  /2 . 59) 

. 74  cm  5. 67  cm 

.9  cm  12.8  cm 

0.78 

•  ((2.59*0. 165) /2 .59) 

.58  cm  5.43  cm 

. 6  cm  12.3  cm 

.79  1.63 

•  ( (2.59-0.418) /2. 59) 

.  1 0  cm  4 . 72  cm 

.  5  cm  1 0 . 7  cm 


3 

) 

5 

2 

0 

) 

5 

2 

0 

) 

5 

1 


%  Error 


2.  61 


4.67 
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Tree  3749 


DB H=7 • 3  cm 
Bark=0.53  cm 
Age=49  yr 
Core=3. 14  cm 


Age  at  Breast  Height=20  yr 
Core  Dif f erential=1 . 456  cm 

Density=1 8 • 66  trees/100  sg.  m 


TC=3 • 14-1.465 
=1.68  cm 


2.18  cm 


1.18  cm 


IR  =  (-0.  16  0+  (0.456*7.  3)  )  •  (  (  1 . 68-0.  1 1  0)  /1 . 68) 

=2.95  cm  3.00  cm  2.86  cm 


DBH  =6.8  cm 


6.9  cm 


6.7  cm 


%  Error  1.45  1.49 

IR  =(-0.  16  0+  (0.456*7.3)  )•(  (1 .68-0.  187)  /I  .68) 

=2.81  cm  2.89  cm  2.66  cm 

DBH  =6.5  cm  6.7  cm  6.2  cm 

%  Error  2.99  4.84 

IR  =(-0. 160+ (0. 456*7.3) ) • ( (1 .68-0.253) /I .68) 

=2.68  cm  2.79  cm  2.48  cm 

DBH  =6.3  cm  6.5  cm  5.8  cm 


%  Error 


3.08 


8.62 


. 


'5  '  . 


■ 

*  i  .  ;  ♦ 


.  ■  • 


126 


Appendix  2.  Multivariate  Statistical  Techniques 

Multiple  Regression 

Multiple  regression  is  a  statistical  technique  designed 
to  examine  the  extent  of  mathematical  dependence  of  a 
particular  variable  (dependent  or  criterion  variable;  Y)  on 
a  set  of  associated  variables  (independent  or  predictor 
variables;  X,  . ..Xn).  The  outcome  of  an  analysis  of  this 
type  is  an  equation  of  the  form 

a  +  fo,  X,  +  b2  X2  "S’,  •  •  bn  Xn 

where  Y  is  the  predicted  value  of  the  dependent  variable 
under  the  specified  conditions  of  the  independent  variables, 
X(  through  Xn  and  a  and  b  •  .  .bn  are  constants  derived  by  the 
procedure  of  least  squares.  This  equation  is  asymmetrical 
in  the  sense  that  although  changes  in  the  magnitude  of  the 
independent  variables  can  be  expected  to  elicit  changes  in 
the  dependent  variable,  the  reverse  is  not  true. 

The  data  must  conform  with  a  series  of  assumptions  to 
be  suitable  for  multiple  regression  analysis  (Poole  1974). 
The  dependent  variable  must  be  normally  distributed  and  the 
relationships  between  it  and  the  independent  variables  must 
be  linear  and  additive.  Although  multiple  regression  is 
quite  robust,  extreme  correlations  between  the  independent 
variables  can  affect  the  results.  Finally,  the  differences 
between  the  observed  and  predicted  values  of  the  dependent 
variables  (residuals)  must  be  normally  distributed  with  a 
mean  of  zero  and  a  constant  variance. 
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Since  regression  analysis  is  associated  with  an 
esoterica  that  often  obscures  its  usefulness,  a  list  of  more 
commonly  encountered  terms  is  provided.  The  general  form  of 
the  multiple  regression  equation  given  above  contained  two 
types  of  constants: 

a  -  represents  the  predicted  value  of  the  dependent 
variable  when  the  values  of  all  the  independent 
variables  are  zero  (Poole  1974), 
b,  .  •  .  b„  -  the  partial  regression  coefficients  represent 
the  magnitude  and  direction  of  change  in  the  dependent 
variable  when  the  associated  independent  variable  is 
increased  by  one  unit  and  there  is  no  change  in  the 
size  of  the  other  independent  variables.  If  these 
coefficients  are  calculated  from  standardized  data 
(transformed  to  have  a  mean  of  zero  and  a  standard 
deviation  of  one,  the  unit  of  measure  being  one 
standard  deviation)  they  are  known  as  beta  weights 
(Draper  and  Smith  1966).  Confidence  in  a  partial 
regression  coefficient  is  expressed  in  its  standard 
error  which  describes  the  variability  involved  in 
sampling  a  particular  population  (Kim  and  Kohout  1975) . 
The  usefulness  of  a  calculated  equation  can  be  measured  in 
several  ways: 

multiple  correlation  coefficient  (R)  -  describes  the 

mutual  relationship  between  the  observed  value  of  the 
dependent  variable  and  the  value  predicted  from  the 
prevailing  values  of  the  independent  variables 
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(Tatsuoka  1969) .  The  square  of  this  value,  R2 
(coefficient  of  multiple  determination)  ,  represents  the 
proportion  of  the  total  variation  of  the  dependent 
variable  explained  by  the  calculated  regression 
equation  (Draper  and  Smith  1966). 
standard  error  of  the  estimate  -  the  dimension  of  one 
standard  deviation  for  observed  values  of  the  dependent 
variable  about  the  regression  equation  (Draper  and 
Smith  1 966) . 

overall  F-test  -  a  significance  test  that  assesses  the 
effectiveness  of  the  regression  equation  to  explain 
more  variation  in  the  dependent  variable  than  would  be 
expected  if  the  equation  were  drawn  at  random  (Draper 
and  Smith  1966) . 

Although  a  particular  equation  offers  a  statistically 
significant  explanation  of  the  variation  in  the  dependent 
variable,  the  importance  of  individual  independent  variables 
could  be  negligible.  Examination  of  each  partial  regression 
coefficient  by  a  partial  F-test  demonstrates  the  ability  of 
the  corresponding  independent  variable  to  explain  variation 
in  the  dependent  variable  left  unexplained  by  other 
variables  in  the  equation. 

Before  attempting  to  calculate  any  regression  equations 
the  data  were  examined  to  ensure  that  all  assumptions  were 
satisfied.  If,  after  conducting  a  Kolmogorov-Smirnov 
goodness  of  fit  test,  the  dependent  variable  was  found  not 
to  be  normally  distributed,  a  suitable  transformation  was 


.  ( 

- 

.  I 


- 

■ 


- 

XJ 


. 

•, 7  56  ad* 


129 


applied.  Simple  correlation  coefficients  were  used  as  an 
index  of  a  linear  relationship  between  the  dependent 
variable  and  five  forms  of  each  independent  variable 
(untransformed,  square,  square  root,  log,  inverse) •  The 
combination  with  the  largest  coefficient  was  selected  for 
use  in  the  analysis.  Following  selection  of  the  "best” 
regression  equation,  the  predicted  values  of  the  dependent 
variable  were  plotted  against  the  actual  values  and  these 
plots  were  visually  examined  for  patterns  that  could  expose 
violations  of  the  assumptions  regarding  the  distribution  of 
residuals  and  the  additivity  requirement. 

A  forward  selection  procedure  was  used  to  derive 
regression  equations  that  contained  only  those  independent 
variables  which  were  useful  in  explaining  variation  in  the 
dependent  variable.  Independent  variables  were  selected  to 
enter  the  equation  one  at  a  time  on  the  basis  of  their 
partial  correlation  coefficient.  This  coefficient  describes 
the  association  between  the  dependent  variable  and  a 
specific  independent  variable  when  the  effects  of  the 
independent  variables  already  in  the  equation  are  held 
constant.  The  process  of  adding  independent  variables  with 
the  highest  partial  coefficients  continued  until  none  of  the 
remaining  potential  independent  variables  could  explain  a 
statistically  significant  amount  of  the  remaining 
unexplained  variation  as  determined  by  partial  F-tests  (Kim 
and  Kohout  1975).  Draper  and  Smith  (1966)  however  suggest 
that  the  forward  selection  technique  is  not  completely 
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adequate  since  a  variable  which  entered  the  equation  at  an 
early  stage  may  become  superfluous  because  of  its 
relationship  with  variables  that  entered  the  equation  later. 
For  this  reason,  each  of  the  partial  F-tests  for  all 
variables  in  the  final  equation  were  examined  for 
significance.  If  the  contribution  of  any  variable  had 
fallen  below  the  acceptance  level  (p  <  0.05)  it  was  removed 
from  the  equation  and  the  selection  process  was  resumed. 

I4dependent  variables  in  the  equation  may  be  highly 
correlated  with  unmeasured  variables  which  are  in  fact 
responsible  for  variation  in  the  dependent  variables 
(Box  1966) .  Because  the  occurrence  of  these  "latent” 
variables  is  to  be  expected  in  the  analysis  of  biological 
systems,  Austin  (1971)  has  stated  that  the  success  of  a 
multiple  regression  equation  is  measured  not  only  by  the 
size  of  the  coefficient  of  multiple  determination,  but  also 
by  the  ability  to  provide  a  biological  interpretation  for 
the  terms  in  the  equation. 
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Discriminant  Analysis  (see  Tatsuoka  1970  and  Morrison  1974) 
The  objective  of  discriminant  analysis  is  to  devise  a 
classification  scheme  that  will  best  distinguish  between 
several  exclusive  groups.  This  technique  involves  the 
derivation  of  a  series  of  equations  that  can  be  used  to 
calculate  a  score  (Z.  )  for  each  individual  or  object  based 
on  the  values  of  a  set  of  independent  variables  (X,.  ...Xni  ) 
associated  with  it  (Morrison  1974).  These  equations  take 
the  form 

Z,  =a,  +b,.  X„  +,..bni  Xni 

in  which  a  and  b,  .  „.b  are  constants  analogous  to  the 

i  1 1  n  i  3 

various  regression  coefficients.  Unlike  regression  analysis 
however,  the  value  of  the  discriminant  score  does  not  bear 
an  exact,  intuitive  connotation;  it  is  merely  a 
classification  scale.  The  probability  of  an  individual  with 
a  given  discriminant  score  belonging  to  one  of  the  exclusive 
groups  varies  along  this  classification  scale.  By  comparing 
each  individual  score  with  a  critical  value  designating  the 
probability  boundary  between  two  groups,  each  individual  can 
be  allotted  to  a  group.  Because  discriminant  functions  are 
linear  they  may  not  adequately  classify  three  or  more  groups 
and  additional  equations  may  be  necessary.  Discriminant 
analysis  was  used  to  distinguish  cnly  two  groups  in  this 
study  (used  and  unused  trees)  so  that  the  complications 
involved  in  more  detailed  analysis  are  not  considered  here 
(consult  Tatsuoka  1970)  . 
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Discriminant  analysis  can  be  used  in  two  ways; 
classification  and  interpretation.  Classification  involves 
deriving  a  series  of  discriminant  functions  from  known 
individuals  and  using  them  to  categorize  unknown  individuals 
on  the  basis  of  associated  independent  variables. 
Interpretation,  the  use  to  which  this  technique  was  put  in 
this  study,  involves  determining  which  of  the  independent 
variables  are  most  effective  in  grouping  individuals.  The 
ability,  of  an  independent  variable  to  increase  the  power  of 
a  discriminant  function  is  inversely  measured  by  Wilks* 
lambda,  which  can  be  transformed  into  a  chi-square  statistic 
for  significance  testing  (Klecka  1975).  Using  a  stepwise 
procedure,  those  independent  variables  that  make  a 
significant  contribution  to  the  discriminant  function  can  be 
identified . 
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Appendix  3.  Detailed  Multiple  Regression  Statistics 

The  statistics  that  follow  are  more  extensive 
descriptions  of  the  relationships  between  variables  outlined 
in  Tables  5,  6,  7,  17  and  23.  The  name  of  the  dependent 
variable  is  found  in  the  upper  left  corner  of  each  block  of 
statistics,  followed  by  the  test  statistic  (Z)  and 
probability  level  for  the  Kolmogorov-Smirnov  test  for 
normality.  The  first  column  of  numbers  below  this 

j 

information  includes  the  statistics  for  the  overall 
regression  equation:  R2  -  the  proportion  of  variation 
explained  by  the  entire  equation:  se  -  the  standard  error  of 
the  estimate;  F  -  the  overall  F-test;  and  cons.  -  the  value 
of  the  Y-intercept.  Columns  to  the  right  contain  statistics 
particular  to  the  respective  independent  variables:  se  -  the 
standard  error  of  the  regression  coefficient;  F  -  the 
partial  F-test;  cons.  -  the  regression  coefficient;  and 
beta  -  the  standardized  regression  coefficient.  Independent 
variables  are  listed  from  left  to  right  in  the  order  in 
which  they  entered  the  regression  equation.  The  R2  value 
associated  with  each  independent  variable  is  a  measure  of 
its  contribution  in  explaining  the  variation  that  remained 
unexplained  by  the  independent  variables  to  its  left.  The 
numbers  in  braces  refer  to  the  type  of  transformation 
applied  to  a  variable:  {1}  squared;  {2}  square  root; 

{3}  log;  {4}  inverse;  {5}  arcsin. 
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Zoratti 

Windward 

Fir  Tables  5,  6  and  7 

DBH  {3} 

Z=0. 481 

p>0. 9 

Age  {4} 

n 

25 

R 

0.  740 

R2 

0.548 

0.548 

se 

0.284 

7.215 

F 

27.874 

27.874 

P 

p<0. 001 

p<0 . 001 

cons. 

1. 591 

-38.094 

beta 

-0. 740 

Height 

{3}  Z=  0.  460 

p>0  •  9 

Age  {4} 

n 

24 

R 

0.632 

R2 

0.399 

0.399 

se 

0.  175 

4.  463 

F 

14.600 

14.600 

P 

p<0. 001 

pCO.QOI 

cons. 

0.898 

-17.054 

beta 

-0.632 

Cl 

Z=0.551  p>0„9 

Age  {4} 

Slope  {4} 

n 

24 

R 

0.757 

R2 

0.573 

0.485 

0.088 

se 

0.670 

17. 131 

0.971 

F 

14.  088 

22.238 

4.335 

P 

p<0.001 

p<0 . 001 

0. 05>p>0.0 

cons. 

3.393 

-80.734 

2.022 

beta 

-0.674 

0.298 
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Skelton  Windward  Fir  Tables  5,  6  and  7 


DBH  {2} 

Z=0. 359 

p>0.  9 

n 

54 

Age  {4} 

Exposure  {1} 

Density  {4} 

R 

0.  736 

R2 

0. 542 

0.217 

0.251 

0.074 

se 

0.714 

24. 733 

0.00001 

2.221 

F 

19.695 

35.592 

31 .509 

8.084 

P 

cons. 

p<0. 001 

2. 769 

p<0.001 

-147.556 

*— 

o  o 
o  o 

•  o 

o  o 

V  • 

0.0 

Q,Q1>p>0.005 

6.316 

beta 

-0.583 

0.553 

0.276 

Height 

{2}  Z=0  • 

786  p>0. 5 

Exposure  {1} 

Age  {2} 

n 

49 

R 

0.552 

R2 

0.304 

0.125 

0.179 

se 

0.  493 

0.000 

0.072 

F 

10. 066 

13.452 

1 1.854 

P 

p<0.001 

p<0.001  0. 

005>p>0 .00 1 

cons. 

-0.729 

0.00002 

0.249 

beta 

0.467 

0.438 

Cl  {3} 

Z=0 • 638 

p>0 .75 

Density  {1} 

Age  {4} 

n 

49 

R 

0.510 

R2 

0.260 

0.149 

0.111 

se 

0.095 

0.0002 

3.641 

F 

8.  067 

10.732 

6.901 

p  0 • 005>p>0. 001  0 

. 005>p>0 .001 

0.025>p>0.01 

cons. 

0.  580 

-0. 0006 

-9.563 

beta 

-0.417 

-0.335 

:  . 
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Tables  5r  6  and  7 


Zoratti  Leeward  Pure  Fir 


DBH  {2}  Z=0 . 589  p>0.75 


n 

104 

Age  {3} 

R 

0.  796 

R2 

0.633 

0.499 

se 

0.  723 

0.502 

F 

87.071 

102.  123 

P 

p<0. 001 

p<0.001 

cons. 

-3.974 

5.072 

beta 

0.624 

Density 


0.133 

0.014 

36.701 

p<0.001 

-0.084 

-0.374 


Height 

Z=1 .309 

p>0. 05 

Age  {3} 

n 

88 

R 

0.  772 

R2 

0.596 

0.522 

se 

2.  849 

2.054 

F 

62.668 

89.919 

P 

p<0.001 

p<0  •  0 

cons. 

-21.228 

19.474 

beta 

0.667 

Density 


0.074 

0.060 

15.501 

p<0.001 

-0.238 

-0.277 


Cl  {2}  Z=0 .6  06  p>0. 7 5 


n 

88 

R 

0.515 

R2 

0.265 

se 

0.  197 

F 

15.326 

P 

p<0. 00 

cons. 

beta 

1.813 

Age  {4} 


Density  {3} 


0.171 

2.787 

13.957 

p<0.001 

-10.411 

-0.354 


0.094 

0.106 

10.823 

0 .005>p>0.001 
-0.347 
-0.312 
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Skelton  Leeward  Pure  Fir  Tables  5,  6  and  7 

DBH  {2}  Z=0 . 566  p>0.9 


Age  {3} 

Density  {4} 

n 

247 

R 

0.  755 

R2 

0.570 

0.487 

0.083 

se 

0.7  33 

0.370 

1.098 

F 

161.590 

277.211 

47.062 

P 

p,<0. 001 

p<0.001 

p<0 .001 

cons. 

”7.525 

6.168 

7.534 

beta 

0.699 

0.288 

Height 

Z  =  1 . 68  0 

p=0.007  |not 

normally  distr 

Age  {3} 

Exposure  {4} 

n 

219 

R 

0.698 

R2 

0.487 

0.476 

0.012 

se 

3.025 

1.487 

735.078 

F 

102. 705 

204.909 

4.913 

P 

p<0. 001 

p<0.001 

0. 05>p>0.025 

cons. 

-21.765 

21.289 

”1629.386 

beta 

0.701 

”0.109 

Cl  {3} 

Z=0 . 8  45 

p>0. 2  5 

Age  {4} 

Density  £3} 

n 

219 

R 

0.656 

R2 

0.431 

0.293 

0.138 

se 

0.116 

1.015 

0.040 

F 

81. 723 

120.538 

52.340 

P 

p<0. 001 

p<0.001 

p<0 .001 

cons. 

0.679 

-11.143 

”0.288 

beta 

”0.565 

-0.372 
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Zoratti 

Leeward 

Mixed  Fir 

DBH  {3} 

Z=0. 929  p>0 .25 

n 

28 

Age  {2} 

R 

0.823 

R2 

0.677 

0.677 

se 

0.271 

0.032 

F 

54.  537 

54.537 

p 

cons. 

p<0.001 

-0.652 

p<0.001 

0.239 

beta 

0.832 

Height 

{3}  Z=0. 660 

p>0. 75 

Age  {3} 

n 

25 

R 

0.792 

R2 

0.628 

0.628 

se 

0.  166 

0.  140 

F 

38.  788 

38.788 

P 

p<0. 001 

p<0.001 

cons. 

-0.695 

0.874 

beta 

0.792 

Cl  {2} 

Z=0 . 672 

p>0. 7  5 

Age 

n 

25 

R 

0.829 

R2 

0.687 

0.687 

se 

0.  182 

0.002 

F 

50.395 

50.395 

P 

p<0. 001 

p<0.001 

cons. 

0.  756 

0.013 

beta 

0.829 

Tables  5, 


6  and  7 


• 

. 

. 

. 

. 

M  . 

' 

| 

. 

Zoratti  Windward  Pine  Tables  5,  6  and  7 


DBH  {2} 

Z=0. 469 

p>0 . 9 

-  w  — - 

n 

94 

Age  {2} 

Slope  {4} 

R 

0.685 

R2 

0.469 

0.425 

0.044 

se 

0.646 

0.047 

0.231 

F 

40.211 

75.274 

7.534 

P 

cons. 

p,<0.  00  1 

-0. 71 1 

p<0.001 

0.412 

0. 0 1>p>0 .005 
0.634 

beta 

0.664 

0.210 

Height 

{3}  Z  =  0 . 80  7 

p>0 . 5 

Age  {3} 

Slope  {4} 

n 

83 

R 

0.597 

R2 

0.356 

0.292 

0.064 

se 

0.  125 

0.  101 

0.047 

f 

22. 137 

39.408 

7.984 

P 

p<0.001 

p<0.001 

0. 01>p>0.005 

cons. 

-0. 728 

0.636 

0.133 

beta 

0.566 

0.255 

Cl 

Z=  0. 960  p>0 . 25 

Age  {2} 

Slope  £1} 

n 

83 

R 

0.678 

R2 

0.459 

0.  420 

0.039 

se 

0.  730 

0.056 

0.0003 

F 

33.964 

65. 556 

5.836 

P 

p<0.00  1 

p<0.001 

0. 05>p>0. 025 

cons. 

-0.842 

0.450 

-0.0008 

beta 

0.670 

-0.200 

. 

. 

. 

. 
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Skelton  Windward  Pine  Tables  5, 


DBH 

Z=1 .044 

►a 

V 

o 

• 

Age 

W 

n 

19 

R 

0. 

652 

R2 

0. 

425 

0. 

425 

se 

0. 

043 

1. 

772 

F 

12, 

585 

12. 

585 

p  0 . 

005>p>0. 001  0. 

005>p>0.001 

cons 

0. 

010 

6  • 

287 

beta 

0. 

652 

Height 

{3}  Z=0. 

803  p>0.5 

Exposure  {1} 

n 

12 

R 

0.  720 

R2 

0.519 

0.519 

se 

0.084 

0.000 

F 

10.780 

10.780 

p  0.025>p>0.01 

0. 025>p>0. 0 1 

cons. 

0.327 

0.00004 

beta 

0.720 

Cl 

Z=0 .4 

98  p>0  •  9 

Density  {4} 

n 

12 

R 

0. 

585 

R2 

0. 

343 

0.343 

se 

0. 

439 

8.008 

F 

5. 

211 

5.211 

p  0. 

05>p> 

0.025  0. 

05>p>0. 025 

cons. 

4. 

554 

-18.280 

beta 

-0. 585 

6  and  7 


. 

. 
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Zoratti  Leeward  Mixed  Pine  Tables  5,  6  and  7 

DBH  Z  =  0. 53  0  p>0 . 9 


n 

31 

Age 

R 

0.678 

R2 

0. 459 

0.459 

se 

5.  296 

0.048 

F 

24.623 

24.623 

P 

p<0.001 

p<0.001 

cons. 

0.472 

0.238 

beta 

0.678 

Height 

{2}  Z=0 . 711 

p>0 . 5 

n 

26 

Age 

Slope  {1} 

R 

0.  742 

R2 

0.511 

0.443 

0.108 

se 

0.359 

0.052 

0.0004 

F 

14.  099 

27.831 

5.538 

P 

cons. 

p<0.001 

0.  932 

p<0.001 

0.274 

0 • 05>p>0. 025 
-0.001 

beta 

0.878 

-0.391 

Cl 

Z=0 . 4  65  p>0 . 9 

Age 

Slope  {1} 

Density  {4} 

n 

26 

R 

0.804 

R2 

0.647 

0.427 

0.110 

0. 110 

se 

0.585 

0.007 

0.001 

2.238 

F 

13. 451 

36.  056 

9.410 

6.888 

P 

p<0.001 

p<0.001 

0 • 005>p>0 • 001 

0 • 005>p>0 .001 

cons. 

0.  877 

0.041 

-0.002 

5.874 

beta 

0.951 

-0.490 

0.338 

. 


. 


. 


. 

. 

, 

. 

. 
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Zoratti  Windward 


Table  17 


Proportion  of  Trees  Used  Z=0.297  p>0.9 


n  9 

R  0.744 

R2  0.553 

se  0.075 

F  8.656 

p  0. 025>p>0. 01 
cons.  0.043 

beta 


Density  {4} 


0.553 
0.  744 
8.656 

0.0  25>p>0. 0 1 
2.  189 
0. 744 


Skelton  Windward 

Proportion  of  Trees  Used  Z=0.618  p>0.75 

Not  significant  F=4.250  p>0.05 


. 


. 
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Zoratti  Leeward  Pure  Table  17 


Proportion  of  Trees 

Used  {5} 

2=0.264  p>0 . 9 

n 

12 

Density  {1} 

Slope  {4} 

R 

0.706 

R2 

0.499 

0.  185 

0.314 

se 

9.011 

12.950 

549.894 

F 

4.482 

8.222 

5.633 

p  0 . 

cons. 

05>p>0. 025  0. 

129.933 

025>p>0. 01 
-37. 132 

0. 05>p>0.025 
-1305.128 

beta 

-0.833 

-0.690 

Skelton 

Leeward  Pure 

Table  17 

Proportion  of  Trees 

Used  Z=0 . 897 

p>0 . 3 

Density  {3} 

Slope  {3} 

Exposure  £1} 

n 

29 

R 

0.  84  7 

R2 

0.717 

0.407 

0.230 

0.080 

se 

0.  136 

0.125 

0.173 

0.000 

F 

21.117 

20.533 

12.336 

7.099 

P 

p<0.001 

p<0.001  0. 

005>p>0 .001 

0 • 025>p>0 . 0 1 

cons. 

-0.286 

-0.564 

0.607 

0.000003 

beta 

-0.500 

0.396 

0.307 

Zoratti  Leeward  Mixed  Table  17 

Proportion  of  Trees  Used  Z=0.457  p>0.9 


n 

5 

Proportion  of  Douglas 

R 

0.958 

R2 

0.917 

0.917 

se 

0.077 

0.145 

F 

33. 162 

33. 162 

p  0  • 

025>p>0. 01 

0.0  25>p>0 . 0 1 

cons 

.  -0.001 

0.835 

beta 

0.  958 

. 
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Zoratti 

Combined 

Table  17 

Proportion  of  Trees  Used  Z=0.470 

p>0 . 9 

#f ir/#total 

X  DBH  {2} 

n 

27 

R 

0.  844 

R2 

0.  713 

0.652 

0.061 

se 

0.  146 

0.  120 

0.060 

F 

29.782 

5,737 

5.073 

P 

p<0. 001 

0. 05>p>0,025  0. 

05>p>0. 025 

cons. 

-0.  193 

0.287 

0.135 

beta 

0.456 

0.429 

Skelton 

Combined 

Table  17 

Proportion  of  Trees 

Used  Z  =  1 . 1 7  0 

p>0. 1 

Exposure  £1} 

Density  £4} 

Slope  £4} 

n 

40 

R 

0.792 

R2 

0.628 

0.348 

0.237 

0.043 

se 

0.151 

0.000 

0.480 

0.916 

F 

20.238 

17.387 

23.422 

4.196 

P 

p<0.001 

p<0.001 

p<0.001 

0.05>p>0. 025 

cons. 

-0.256 

0.000003 

2.322 

-1.875 

beta 

0.479 

0.492 

-0.235 

. 

. 

. 

. 

r  •• 

.  - 

. 

Zoratti 


Table  23 


X  DBH 

Z=0. 59 

n 

17 

R 

0.831 

R2 

0.690 

se 

3.089 

F 

33.337 

P 

p<0.00 

cons. 

30.503 

beta 

X  Height  {1} 

p>0 .75 

Use  Age  {3} 


0.690 

2.363 

33.337 

p<0.001 
-  13. 644 
-0.831 


Z=0. 6  19  p>0. 75 


Not  significant  F=3.815  p>0.05 


X  Cl 

{1}  Z=0 .758 

p>0.  5 

Use  Age  {4} 

n 

17 

R 

0.  841 

R2 

0.707 

0.707 

se 

0.  934 

1.035 

F 

36.  125 

36.  125 

P 

p<0. 001 

p<0. 001 

cons. 

2.  620 

6.222 

beta 

0.941 
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Skelton  Table  23 

X  DBH  {1}  2=0 

.431 

p>  0 . 9 

Not 

significant 

F= 

4.  277  p>0  • 

X  Height  {1} 

Z=0 . 6 

02  p>0. 75 

Not 

significant 

F= 

1. 446  p>0.  1 

X  Cl 

{4}  Z=0 . 

593 

p>0 . 8  74 

Use  Age  {4} 

n 

16 

R 

0.526 

R2 

0.277 

0.277 

se 

0.050 

0.0001 

F 

5.366 

5.366 

P  o 

•  05>p>0. 025 

0. 

05>p>0. 025 

cons 

.  0.409 

0.0003 

beta 

0.526 

. 


